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Significant  results  of  the  second  year  of  this  study  Include  the 
following.  A  commercial  foward-angle  laser  light-scattering  Instrument  has 
been  modified  to  measure  drop  sizes  in  evaporating  and  burning  sprays.  Using 
this  apparatus  it  has  been  shown  that  emulsified  and  neat  fuels  initially 
atomize  to  about  the  same  drop  sizes,  but  at  elevated  temperatures  and 
pressures,  both  Jet-A  and  hexadecane  emulsions  evaporate  in  such  a  way  so  as 
to  produce  drops  of  significantly  smaller  Sauter  mean  diameter  (SMD)  than 
their  neat  fuel  counterparts  at  the  same  distance  from  the  nozzle.  The 
difference  between  the  emulsified  and  neat  fuel  drop  sizes  only  appear  at 
elevated  pressures  (greater  than  about  3-1/2  atm)  for  the  fuels  and 
temperatures  examined.  These  results  are  consistent  with  those  that  might  be 
expected  if  microexplosions  were  occurring  in  the  emulsified  fuels. 


A  computer  model  has  been  developed  based  on  single  drop  heat-up  and 
evaporation  theory.  Sixteen  to  twenty-one  size  classes  are  used  to 
characterize  the  spray  and  the  complete  heat-up,  evaporation,  and  aerodynamic 
interactions  are  included  in  the  model.  The  model  is  Incomplete  in  that  the 
lnhomogene itles  of  the  spray  along  the  line-of-sight  have  not  yet  been 
Included.  However,  the  steady  state  drop  temperature  has  been  computed  and 
compared  with  the  spontaneous  nucleatlon  temperature  of  water  for  a  number  of 
conditions.  These  calculations  suggest  that  the  steady  state  drop  surface 
temperatures  at  conditions  where  differences  are  seen  in  the  evaporation 
characteristics  of  the  emulsified  and  neat  fuels  are  somewhat  below 
(approximately  30K)  the  kinetic  superheat  limit  for  water.  Spontaneous 
nucleatlon  is  expected  to  occur  somewhat  below  the  superheat  limit. 
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MECHANISMS  OF  SMOKE  REDUCT  101  IN  THE  HIGH  PRESSURE 
COMBUSTION  OR  EMULSIFIED  FUELS; 

VOL.  II:  EXPERIMENTAL  AND  THEORETICAL  STUDY  OF  EVAPORATING 
EMULSIFIED  AND  NEAT  FUEL  SPRAYS 


INTRODUCTION  AND  BACKGROUND 

High  quality  petroleum  derived  fuels  are  not  as  widely  available  as 
they  were  In  the  past.  Future  fuels  derived  from  petroleum  and  non¬ 
petroleum  sources  such  as  coal  and  shale-oil  will  tend  to  be  lower  In 
hydrogen  content  and  more  viscous  with  more  residual  components.  These 
fuels  will  be  more  difficult  to  atomize  effectively  and  will  tend  to 
produce  more  soot  during  combustion.  The  purpose  of  this  contract  Is  to 
conduct  Fundamental  research  which  addresses  two  separate  aspects  of  the 
problems  associated  with  atomization  and  burning  of  lower  quality  fuels. 
The  first  aspect  Is  the  experimental  study  of  atomization  and  evaporation 
of  water /oil  emulsified  fuel  sprays  In  a  real  combustion  environment, 
with  a  specific  effort  to  verify  the  proposed  "mlcroexploslon"  phenomena. 
The  second  goat  Is  to  examine  the  atomization  and  evaporation  of  "out-of- 
spec”  fuels  which  are  too  viscous  or  non-volatile  to  meet  the  standard 
fuel  specifications  for  gas  turbine  combustion.  Both  of  these  goals 
required  the  development  of  experimental  techniques  not  previously 
available  to  study  drop-size  distributions  of  fuel  sprays  In  high 
pressure  evaporating  and  combusting  systems.  This  report  covers  progress 
made  during  the  second  year  of  a  three  year  effort  to  study  atomization 
in  combustors. 

The  first  aspect  of  the  atomization  problem  deals  with  the  study  of 
emulsified  fuels.  It  has  been  shown  (Refs.  1  and  2)  that  the  Increased 
soot  formation  resulting  from  combustion  of  low-hydrogen-content  fuels 


can  be  somewhat  offset  by  mixing  water  or  alcohol  with  the  fuel.  Several 
theories  have  been  offered  to  explain  this  soot-reductlon  phenomena. 
These  Include:  (a)  better  atomization  of  the  emulsified  fuel  due  to 
"microexploslons"  caused  by  superheating  and  violent  boiling  of  Interior 
water  micro-droplets;  (b)  a  change  in  the  chemical  kinetics  leading  to 
soot  production  due  to  the  Increased  OH  or  the  increase  in  H/C  ratio;  (c) 
the  reduction  of  liquid-phase  pyrolysis  processes  due  to  lower  drop 
temperatures  resulting  from  the  lower  boiling  point  of  water;  or  (d)  the 
reduction  of  gas-phase  pyrolysis  processes  due  to  the  lower  flame 
temperature  from  either  strictly  dilution  effects  or  other  chemical 
processes.  One  of  the  goals  of  this  contract  was  to  develop  techniques 
to  verify  the  existence  of  microexploslons  in  real  combustion  systems. 
Microexploslons  have  been  observed  for  large  emulsified  drops  burning  in 
quiescent  atmospheres,  but  no  direct  measurements  have  been  successfully 
performed  on  realistic  sprays  in  turbulent  air.  Some  of  the  observed 
microexplosion  phenomena  have  been  criticized  as  unrealistic  for  several 
reasons:  nucleatlon  Introduced  by  the  supporting  wire;  the  large  size  of 
the  drops  (when  compared  with  real  spray  size  distributions);  and  the 
lack  of  high  aerodynamic  shear  forces  which  promote  internal  mixing  and 
reduce  the  possibilities  for  microexploslons  and  also  extend  the  heat-up 
period  relative  to  the  total  lifetime  of  the  drop.  The  purpose  of  the 
first  part  of  this  program  was  to  develop  methods  to  examine  the  behavior 
of  sprays  of  emulsified  fuels  in  a  more  realistic  combustion  environment. 

The  second  part  of  this  program  involves  a  study  of  the  atomization 
and  evaporation  of  alternative  fuels,  which  generally  will  be  more 
difficult  to  atomize  and  evaporate  when  compared  to  present-day 
petroleum-based  fuels.  Some  correlations  exist  for  predicting  the  change 
in  spray  quality  for  small  variations  in  fuel  properties,  and  some  models 
have  been  developed  for  fuel  spray  evaporation  in  combustors,  but 
experimental  data  are  lacking  to  verify  these  models  and  extend  them  to 
alternative  fuels.  The  same  experimental  techniques  developed  for  the 
first  part  of  the  program  will  be  used  for  the  second  part. 


EXPERIMENTAL  APPARATUS 


The  combustor  facility  and  test  section  were  described  in  some 
detail  in  the  previous  yearly  progress  report  (Ref.  3).  That  system  was 
used  for  these  tests  except  that  the  disc,  which  provides  recirculation 
for  flame  stabilization,  was  replaced  with  a  fuel  nozzle  alone  for  the 
evaporation  tests,  and  a  window  purge  was  added.  These  modifications,  as 
well  as  a  brief  description  of  the  overall  system  are  discussed  below. 

Unvitiated  compressed  air  is  supplied  to  the  test  cell  by  three 
compressors  and  a  gas-fired  preheater  which  provide  up  to  1.1  kg/ sec  at 
pressures  ranging  from  138  kPa  (1.36  atm)  to  1620  kPa  (16.0  atm)  and 
temperatures  up  to  1090  K. 

The  test  section  consists  of  a  type  316  stainless  steel  pipe, 
16.83  cm  0D,  12.50  cm  ID,  2.17  cm  wall  thickness,  with  quartz  windows  on 
each  side  having  a  clear  aperture  of  9.2  cm  by  6.7  cm.  A  smooth 
transition  section  5.5  cm  long  is  provided  on  each  end  of  the  windows  to 
reduce  recirculation  zones.  Each  window  is  purged  by  air  on  both  top  and 
bottom.  Air  is  supplied  to  each  window  by  a  total  of  10  tubes  of  nominal 
0.635  cm  00  and  then  diffused  and  exhausted  by  a  convergent  section 
parallel  to  the  window  surface. 

For  the  evaporation  tests,  the  flame  stabilizing  disc  has  been 
replaced  by  a  fuel  nozzle  alone  mounted  about  even  with  the  upstream  edge 
of  the  windows.  The  nozzle  is  supported  with  a  12.7  mm  OD  stainless 
steel  tube  which  follows  the  test  section  centerline  for  a  distance  of 
about  250  cm  upstream  from  the  nozzle  and  then  makes  a  right  angle  bend 
and  exits  from  the  test  section.  The  12.7  mm  0D  tube  carries  cooling 
water  and  a  separate  parallel  6.35  mm  0D  stainless  steel  tube  is  the 
water  return  line.  Another  6.35  mm  0D  stainless  steel  tube  is  mounted 
Inside  the  12.7  mm  water  cooling  tube  and  carries  fuel  to  the  nozzle. 
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Fuel  temperature  Is  regulated  by  controlling  the  water  cooling 
temperature. 

The  nozzles  used  were  simplex  pressure  atomizers  producing  a  hollow 
cone  spray.  The  model  type  was  Delavan  Corp.  WDA  45°  hollow  cone  with 
nominal  flow  rates  of  3.79  liters/hour  (1  gallons/hour)  and  11.3  liters/ 
hour  (3  gallons/hour)  at  862  kPa  (125  psi)  differential  pressure  using 
water. 

The  operation  of  the  test  cell  is  monitored  with  a  Hewlett-Packard 
9820  programmable  calculator  which  is  coupled  to  a  50  channel  scanner.  A 
test  report  including  the  average  and  statistical  variation  of  all  the 
important  parameters  is  available  immediately  after  a  test.  The  sensing 
systems  consist  of  strain-gage  pressure  transducers,  thermocouples,  and 
turbine  flowmeters. 

Drop  size  data  were  obtained  with  a  Malvern  Model  2200  Particle 
Sizer  based  on  the  diffraction  angle  produced  by  drops  when  illuminated 
by  a  beam  of  monochromatic,  coherent,  collimated  light  from  a  HeNe  laser. 
The  smaller  drops  diffract  light  at  larger  angles  to  the  optical  axis 
than  the  larger  drops.  Detection  is  accomplished  with  a  30  annular  ring 
set  of  solid  state  detectors.  Detector  outputs  are  multiplexed  and  the 
data  signal-averaged  with  a  Commodore  PET  computer.  A  fairly  complex 
(and  proprietary  to  Malvern)  computer  routine  is  used  to  interpret  the 
light  scattering  pattern  of  the  polydlsperse  drop  systems.  The  drop-size 
data  are  available  as  either  a  set  of  two  parameters  defining  a  Rosln- 
Rammler,  log-normal,  or  normal  distribution,  or  as  a  histogram  of  15  size 
classes  of  drops  without  any  assumption  about  the  shape  of  the 
distribution.  Although  all  of  the  distributions  have  been  tried,  the 
Rosin-Rammler  has  generally  provided  the  best  fit  and  was  used  throughout 
this  program.  A  300  mm  focal  length  f/7.3  lens  was  used  to  collect  the 
scattered  light.  The  laser  beam  diameter  was  9  mm  with  a  Gaussian 
Intensity  distribution  trunc.  ad  at  tv  edge  by  the  9  mm  aperture. 


Several  modifications  were  made  to  the  Malvern  drop  sizer  for  these 
measurements  to  avoid  problems  associated  with  flame  radiation  during  the 
burning  spray  measurements,  and  laser  beam  steering  problems  during  both 
non-burning  evaporation  tests  and  burning  sprays. 

Three  modifications  were  made  to  discriminate  against  flame 
radiation.  First  the  laser  beam  was  chopped  at  667  Hz  and  the 
differential  signal  at  each  detector  ring  was  sampled  for  20  cycles. 
This  necessitated  slowing  the  multiplexing  rate  in  the  electronics  from 
about  1.3  ms/channel  to  30  ms/channel.  It  was  also  necessary  to  sample 
only  the  "clean"  part  of  the  signal  when  the  beam  was  completely  blocked 
or  completely  transmitting  as  the  chopper  blade  edge  produces  a  large 
spurious  signal  when  passing  through  the  laser  beam  due  to  its  own 
diffraction  pattern.  The  second  modification  was  the  use  of  a  3  nm  band 
pass  Interference  filter  centered  at  the  HeNe  laser  wavelength  of  632.8 
nm  in  front  of  the  lens  on  the  drop  sizer.  All  measurements  were  made 
with  the  300  mm  focal  length  lens  which  has  a  maximum  effective 
acceptance  angle  of  about  2.7°  from  the  normal  axis  for  light  reaching 
the  detector.  The  transmission  of  the  filter  changes  less  than  one 
percent  for  rays  at  this  angle  or  less  relative  to  the  normal  axis 
through  the  filter.  However,  if  the  100  mm  or  63  mm  focal  length  lenses 
(the  other  standard  sizes)  were  used,  the  acceptance  angle  Increases  to 
about  8.1°  or  12.8°,  respectively,  and  wider  band  pass  filters  would  be 
necessary.  The  third  modification  Involves  the  use  of  an  aperture  of 
diameter  32  mm  located  about  130  mm  from  the  lens  to  block  radiation  from 
the  flame  which  is  not  within  the  scattering  volume  of  the  laser  beam. 

In  addition  to  problems  with  flame  radiation,  there  are  substantial 
problems  even  in  non-combusting  sprays  in  high  temperature  air  due  to 
steering  of  the  laser  beam  caused  by  refractive  index  changes  in  the 
scattering  volume  resulting  from  thermal  and  compositional  gradients  near 
the  drops.  This  results  in  excess  spurious  signals  on  the  inner  detector 
rings,  getting  worse  at  higher  evaporation  rates  and  finally  limiting,  at 
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some  distance  from  the  nozzle,  the  application  of  this  technique  for  drop 
sizing.  A  technique  has  been  developed  for  blocking  light  on  the  inner 
detector  rings  and  extrapolating  values  from  the  outer  rings.  Radiation 
data  on  up  to  eight  of  the  thirty  rings  is  blocked,  but  even  this  becomes 
insufficient  at  some  point  in  the  evaporation  process. 

The  significant  fuel  property  data  are  given  in  Table  1.  The 
^mulsions  were  first  blended  using  a  homogenizer,  and  were  also 
continuously  recirculated  through  an  ultrasonic  cavity  during  use  to 
maintain  a  consistent  blend. 
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Table  1.  Some  Properties  or  Fuels  Used 


JET  A  (AL-10112F) 

Viscosity  at  40°C  1.72  cS 
Specific  gravity  at  60°F  0.8049 

Ileat  of  Combustion  46.904  MJ/kg 
Hydrogen  Content  14.15Z 
Flash  Point:  355K 
Boiling  Point  Distribution  (GC) 


Wt%  Off 


Deg  K 


0.1 

436 

10 

468 

20 

484 

30 

490 

40 

494 

50 

502 

60 

506 

70 

510 

80 

518 

90 

527 

99 

548 

100 

796 

20Z  H20/2Z  Surf act ant/ 7 8Z  Jet  A  (AL-10112F) 

Viscosity  at  40°C  4.47  cS 

Specific  gravity  at  60°F  0.8618 


Hexadecane  (95Z  purity) 

Viscosity  at  40oc  2.98  cS 
Specific  gravity  at  60°F  0.7774 

Boiling  point  (nominal)  560K 


20Z  H2O/2Z  Surfactant/78%  Hexadecane 


Viscosity  at  40°C  5.55  cS 

Specific  gravity  at  60<>F  0.8249 


SPRAT  TRAJECTORY  AND  EVAPORATION  MODEL 

A  computer  model  based  on  an  array  of  different-sized, 
noninteracting  drops  has  been  developed  to  study  the  heat-up,  final  "wet- 
bulb"  drop  temperature,  evaporation  rates,  trajectories  and  resulting 
size  distribution  parameters  at  locations  downstream  from  the  nozzle. 
The  trajectory  part  of  the  model  should  be  considered  as  a 
straightforward  and  simple  one  based  on  an  average  drop  path  (neglecting 
turbulent  statistical  fluctuations)  in  a  one-dimensional  flow  field. 
This  corresponds  to  the  experimental  situation  being  studied.  The 
thermodynamic  part  of  the  model,  however,  is  very  detailed  and  computes 
the  transient  heat-up  and  final  drop  temperature  in  greater  detail  than 
most  spray  models.  The  environment  surrounding  the  drop  is  calculated 
based  on  the  properties  of  both  the  fuel  vapor  and  air  rather  than  air 
alone  as  in  some  spray  models.  The  convective  effects  on  evaporation  are 
included  and  constantly  updated  as  the  different  size  drops  are 
decelerated  to  the  air-stream  velocity.  The  effects  of  temperature  and 
pressure  on  air  and  fuel  properties  is  Included,  and  the  steady-state 
"wet-bulb"  temperature  is  iteratively  calculated  using  the  equations  for 
the  vapor  pressure  and  temperature-dependent  fuel  and  air  properties. 
This  model  assumes  uniform  drop  temperatures  (infinite  heat  transfer 
rates  within  the  drop)  and  single-component  fuels. 

Some  spray  models  have  predicted  the  change  in  drop-size 
distribution  for  a  quiescent  spray  as  a  function  of  time,  but  results 
from  this  model  show  the  considerable  distinction  between  drop-size 
distributions  as  measured  at  a  given  location  along  the  nozzle  axis,  at 
which  drops  arrive  after  different  transit  times  (depending  on  diameter), 
as  contrasted  with  the  variation  with  time.  For  example,  in  this 
experiment,  the  spray  initially  has  a  higher  axial  velocity  than  the  air 
and  the  small  drops  decelerate  more  rapidly  than  the  large  ones, 
resulting  in  a  higher  concentration  in  the  measurement  volume  for  the 
small  drops  and  a  longer  time-of-arrtval .  These  effects  must  be  included 


in  any  model  which  is  intended  to  be  used  for  comparison  with 
experimentally  measured  changes  in  size  distribution. 

In  developing  a  model  which  can  be  used  to  compare  with  experimental 
results,  it  is  necessary  to  consider  the  interaction  of  the  measuring 
device  with  the  spray.  The  added  weighting  of  the  small  drops  due  to 
their  more  rapid  deceleration  has  been  discussed.  In  addition,  the  small 
drops  are  caught  up  in  the  air  stream  and  tend  to  populate  the  center  of 
the  spray  while  the  large  drops  maintain  their  initial  velocity  longer 
and  tend  to  fill  the  outer  edges  of  the  spray.  This  disparity  is 
magnified  by  the  fact  that  at  the  exit  of  a  hollow-cone  swirl  nozzle  (the 
type  used  in  these  tests)  the  majority  of  the  fluid  consists  of  larger 
drops  in  the  outer  cone  of  the  spray,  while  smaller  drops  fill  in  the 
less  dense  spray  in  the  center  of  the  cone.  Thus,  it  is  necessary  to 
model  these  lnhomogeneitles  in  size  distribution  through  the  spray  and 
also  the  location,  diameter,  and  intensity  distribution  of  the  laser  beam 
relative  to  the  spray.  At  this  time,  the  model  does  account  for  the 
different  trajectories  of  the  large  and  small  drops,  but  assumes  that  all 
the  spray  is  initially  at  a  fixed  cone  angle.  The  drop  size  distribution 
is  calculated  for  a  complete  cross  section  of  the  spray  at  a  given  axial 
location,  without  regard  for  the  laser  beam  sampling  volume.  The  model 
also  assumes  that  spray  breakup  is  complete  and  no  secondary  atomization 
or  drop  agglomeration  occurs  beyond  some  point  from  the  nozzle  exit.  An 
initial  size  distribution  at  that  point  is  an  input  to  the  model.  Thus, 
the  model  is  not  applicable  very  close  to  the  nozzle  (  10mm). 

* 

Results  from  the  model  are  presented  in  the  following  sections, 
while  the  detailed  mathematics  of  the  model  are  presented  in  Appendix  A. 
A  general  flow  chart  of  the  model  is  presented  in  Figure  A-l.  while  a 
listing  of  the  Fortran  source  code  is  given  in  Appendix  R.  Most  of  the 
model  has  been  constructed  from  pieces  of  existing  models.  The 
thermodynamic  part  of  the  model  follows  closely  the  work  of  Chin  and 
Lefebvre  at  Purdue  University,  while  the  aerodynamic  equations  for  drag, 
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velocity,  etc.  are  based  on  a  model  developed  at  the  University  of 
Sheffield  by  Swithenbank,  Boyson,  Ayers  and  others.  (See  Appendix  A  for 
detailed  references.) 

The  capabilities  of  the  model  can  probably  be  best  understood  by 
examining  some  typical  output  data  and  also  comparing  the  model  results 
with  experimental  results.  Some  information  available  from  the  model  is 
shown  in  Table  2  which  represents  a  spray  of  n-hexadecane  into  air  at  644 
K  and  4.42  atm.  The  air  velocity  was  assumed  to  be  a  uniform  one¬ 
dimensional  axial  flow  of  8.S4  m/s,  while  the  initial  fuel  velocity  was 
20.76  m/s  at  an  angle  of  22.5°  to  the  axis  (45°  cone  angle)  and  a  fuel 
temperature  of  311  K.  The  computed  drop  size  data  assuming  a  Rosin- 
Raramler  distribution  are  shown  on  the  top  line  and  Include  the  two  Rosin- 
Rammler  parameters,  x  and  N  which  specify  a  size  and  a  width  for  the 
distribution,  and  the  Sauter  mean  diameter  (SMD)  which  is  an  "average" 
drop  size.  By  definition,  if  all  drops  in  the  spray  were  equal  in  size 
to  the  SMD,  then  the  volume-to-surf ace-area  ratio  of  the  imaginary 
monodlsperse  spray  would  be  equal  to  the  actual  spray.  This  ratio  is 
chosen  because  surface  area  regulates  evaporation  rates  while  the  volume 
indicates  the  quantity  of  fuel.  The  initial  drop  sizes  are  shown  in 
column  l  and  the  initial  relative  population  of  drops  in  each  size  class 
were  computed  to  correspond  exactly  to  a  Rosin-Rammler  distribution 
specified  as  an  input  to  the  program.  (Those  values  were  x  **  57  and  N  * 
1.5  for  the  example  shown.)  Because  of  evaporation  and  the  more  rapid 
deceleration  of  the  small  drops  (resulting  in  an  increased  weighting  for 
the  small  drops),  the  size  distribution  changes  at  positions  downstream 
from  the  nozzle.  The  drop-size  distribution  no  longer  corresponds 
exactly  to  the  Rosin-Rammler  distribution,  and  the  degree  to  which  it 
does  fit  that  distribution  is  given  by  the  correlation  coefficient.  For 
the  cases  examined,  the  correlation  coefficient  has  been  in  excess  of 
0.97,  but  under  some  conditions  it  could  be  worse.  These  distribution 
parameters  are  calculated  from  information  shown  for  individual  drop  size 
classes  in  the  lower  part  of  Table  2. 
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The  initial  drop  sizes  (in/um)  are  shown  in  column  1  while  the 
instantaneous  drop  sizes  are  shown  in  column  2  at  a  distance  from  the 
nozzle  of  26.4  mm  as  shown  in  column  4  and  a  time  (in  seconds)  as 
specified  in  column  3.  When  the  drops  reach  a  diameter  of  0.9  ym  they 
are  assumed  to  be  completely  evaporated  and  the  diameter  is  fixed  at 
0.9  ym  to  keep  some  of  the  computations  from  becoming  undefined.  For 
those  cases,  the  time  and  location  are  frozen  at  the  point  where  the  drop 
disappeared  and  those  drops  are  not  used  in  calculating  the  Rosin-Rammler 
parameters.  Note  that  this  model  calculates  spray  data  at  a  given  axial 
distance  from  the  nozzle  (X  DROP)  and  not  at  a  given  time  from  injection. 
The  smaller  drops  take  longer  to  reach  the  specified  location  due  to 
their  lower  velocity  as  specified  in  column  7  by  U-VEL  in  m/s.  The 
smaller  drops  decelerate  to  the  air  stream  velocity  of  8.54  m/s  more 
rapidly  than  the  large  drops.  The  coordinate  system  is  cylindrical  with 
the  axial  coordinate  being  x  with  velocity  u,  the  corresponding  radial 
components  y  and  v,  and  the  corresponding  angular  components  (assumed 
zero  presently)  z  and  w.  Distances  are  in  meters  and  velocities  in 
meters /second. 

Note  the  difference  in  trajectories  and  velocities  for  the  different 
size  classes.  All  non-evaporated  drops  are  at  the  same  axial  (X  DROP) 
location  but  the  radial  location  (Y  DROP)  of  the  small  drops  is  less  than 
the  large  ones  due  to  their  being  caught  up  in  the  air  stream  sooner  as 
shown  by  the  axial  (U-VEL)  and  radial  (V-VEL)  velocities.  The  wide 
variation  in  Reynolds  numbers  (RE)  is  shown  in  the  last  column.  The 
Importance  of  convective  effects  in  heat  and  mass  transfer  is  obvious 
from  these  values,  and  hence,  the  shortcomings  Inherent  in  quiescent 
spray  models. 


The  instantaneous  drop  temperature  (FTEMP)  is  shown  in  the  next  to 
the  last  column.  The  initial  fuel  temperature  was  311  K  and  the  "wet- 
bulb'*  final  steady-state  drop  temperature  is  520  K.  at  these  conditions. 
The  normal  boiling  point  of  n-hexadecane  is  560  K.  Note  that  the  small 


heat  up  and  evaporate  almost  Instantaneously  while  the  "uniformly  mixed" 
temperature  of  the  largest  drops  is  predicted  to  be  close  to  the  initial 
fuel  temperature.  Due  to  finite  rate  heat  transfer  the  surface 
temperature  of  the  large  drops  would  be  higher,  but  experimental  evidence 
is  lacking  to  predict  temperature  profiles  within  the  drop.  Complex 
models  have  been  developed  to  predict  these  profiles,  but  they  are  too 
cumbersome  for  this  effort  and  generally  lack  verification. 

Note  the  Importance  of  modeling  the  transient  nature  of  the  drop 
before  it  reaches  the  steady-state  temperature.  Of  the  21  drop  size 
classes  shown  in  Table  2,  only  one  nonevaporated  size  class  has  reached 
the  steady-state  temperature.  As  the  pressure  and  convective  effects  are 
Increased,  the  length  of  the  transient  period  relative  to  the  steady- 
state  period  increases  (Ref.  4).  Thus,  models  that  assume  no  mass 
transfer  until  the  drop  reaches  its  steady-state  temperature  are  not 
representative  of  many  realistic  conditions  in  gas  turbines  and  are 
Insufficient  to  compare  with  experimental  results  generated  in  this 
program. 

In  addition  to  the  drop  array  data  just  discussed,  the  model  has 
been  used  to  examine  the  steady-state  "wet-bulb”  temperature  and  the 
quiescent,  evaporation  rate  of  the  fuel,  which  are  both  independent  of 
drop  size,  as  a  function  of  air  pressure  and  temperature.  At  atmospheric 
pressure,  the  steady-state  drop  temperature  is  lower  than  the  normal 
boiling  point  due  to  the  cooling  effect  of  energy  used  in  vaporizing  the 
fuel.  As  the  pressure  Increases,  the  boiling  point  increases  and  the 
steady-state  temperature  Increases.  Obviously  the  fuel  and  air 
properties  are  Important  in  determining  this  temperature.  Steady-state 
temperatures  for  three  normal  paraffins  -  octane  (CgHig), 
dodecane  (Ci2H26)»  and  hexadecane  (C16H34)  -  were  calculated  over  a  range 
of  conditions  with  the  results  as  shown  in  Figures  1  to  3.  At  low  air 
temperatures,  the  drop  temperatures  converge  to  the  air  temperature.  As 
the  air  temperature  increases,  the  Increased  vaporization  cools  the  drop 
below  the  air  temperature,  but  Increased  pressure  markedly  increases  the 
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FIGURE  2.  COMPUTED  STEADY  STATE  DROP  SURFACE  TEMPERATURE  FOR  N-DODECANE 


allowable  drop  temperature.  The  kinetic  limit  of  superheat  for  water 
Indicated  In  the  figures  Is  discussed  In  a  later  section.  The  critical 
pressure  for  hexadecane  is  1421  kPa  (14.02  atm)  and  for  dodecane  1820  kPa 
(17.96  atm),  and  data  are  not  shown  for  pressures  higher  than  the 
critical  pressure. 

Droplet  evaporation  follows  a  D-squared  law  which  relates  the  drop 
diameter  D  to  the  initial  size  D0  and  the  time  t  by  (Ref.  5) 

D2  -  D02  -  Xt  (1) 

where  the  proportionality  constant  X  is  the  evaporation  constant.  The 
mass  transfer  rate  of  fuel  from  the  drop,  mp,  at  quiescent  conditions  is 
related  to  the  evaporation  constant  by 

Ap  -  (jt/4)  XPF  d  (2) 

where  P  p  is  the  fuel  density.  The  evaporation  constant  is  increased  by 
convection,  with  a  correction  factor  given  by  Frossllng  (Ref.  6)  for 
cases  where  heat  transfer  rates  are  controlling  of, 

^Conv.  "  ^Quies.(l  *  0*276  Rep®*^  Prg®*^) 

where  Rep  Is  the  Reynolds  number  of  the  gas  relative  to  the  drop  and  Prg 
is  the  Prandtl  number  of  the  gas.  The  correction  factor  depends  on  drop 
size  and  must  be  calculated  using  the  drop  array  part  of  the  program,  but 
It  Is  instructive  to  examine  the  temperature  and  pressure  dependence  of 
the  quiescent  evaporation  constant  which  is  independent  of  drop  size  and 
velocity.  The  quiescent  evaporation  constants  for  n-octane,  n-dodecane, 
and  n-hexadecane  are  shown  in  Figures  4  to  6.  It  can  be  seen  that  there 
Is  a  strong  temperature  dependence,  as  there  must  be,  and  a  pressure 
dependence  which  is  weak  over  the  conditions  used  in  these  experiments 
but  which  Increases  at  elevated  temperatures. 


FIGURE  4.  COMPUTED  EVAPORATION  CONSTANT  AT  QUIESCENT 
CONDITION  FOR  N-HEXADECANE 


FIGURE  5.  COMPUTED  EVAPORATION  CONSTANT  AT  QUIESCENT 
CONDITION  FOR  N-DODECANE 


P{ ATM ) 
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FIGURE  6.  COMPUTED  EVAPORATION  CONSTANT  AT  QUIESCENT 
CONDITION  FOR  N-OCTANE 


In  summary,  the  computer  model,  while  not  completely  comparable  with 
experimental  results  in  its  present  form,  does  provide  insight  into  much 
of  the  thermodynamics  affecting  drop  evaporation.  It  also  shows 
quantitatively  how  the  smaller  drops  are  caught  up  in  the  air  stream  much 
more  quickly  than  the  larger  ones.  It  also  indicates  the  Importance  of 
modelling  evaporation  during  the  transient  heat-up  phase,  particularly 
for  realistic  sprays  at  high  pressures  and  Reynolds  numbers.  It  shows 
the  effect  of  air  pressure  on  the  steady-state  drop  temperature,  and  how 
increasing  air  pressure  leads  to  increased  drop  temperatures  which 
enhances  the  probability  of  initiating  microexplosions. 
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EXPERIMENTAL  RESULTS  &  DISCUSSION 


Experiments  determining  the  change  of  drop  size  distribution  with 
distance  from  the  nozzle  for  emulsified  and  neat  Jet  A  and  n-hexadecane 
In  high  pres  sure/ temperature  flowing  air  have  been  performed.  Initial 
experiments  were  performed  with  Jet  A  and  showed  that  the  emulsified 
fuels  Initially  atomize  to  about  the  same  size  as  the  neat  fuels  and 
evaporate  In  a  similar  manner  at  low  temperature  and  pressure  conditions. 
However,  at  elevated  temperatures  and  pressures,  the  emulsified  fuels 
evaporated  in  such  a  way  so  as  to  develop  smaller  drop  sizes  than  the 
neat  fuel.  Because  Jet  A  is  a  multi-component  fuel,  as  are  all 
commercial  fuels,  the  evaporation  of  even  neat  Jet  A  requires  several 
Important  assumptions  about  the  distillation  process  to  model.  For  that 
reason,  further  tests  were  conducted  with  a  single  component  fuel,  n- 
hexadecane,  and  the  evaporation  process  for  a  spray  of  that  fuel  has  been 
modelled  and  Is  compared  with  experimental  results. 

The  experimental  results  presented  in  this  section  all  have  the  same 
general  form  consisting  of  graphs  showing  the  average  drop  size, 
represented  by  the  Sauter  mean  diameter  (SMD),  upon  initial  atomization 
close  to  the  nozzle,  and  after  varying  amounts  of  evaporation  moving 
downstream  from  the  nozzle.  In  most  of  the  plots,  the  SMD's  were 
normalized  near  the  nozzle  so  that  the  relative  effects  of  evaporation  on 
the  emulsified  and  neat  (without  water)  fuels  could  be  seen  more  clearly. 
Each  graph  shows  the  change  in  average  drop  size  as  a  function  of 
distance  from  the  nozzle  for  a  series  of  air  temperatures,  all  at  a  fixed 
air  pressure,  or  visa  versa.  Data  for  neat  fuels  are  shown  as  open 
symbols  with  solid  lines  drawn  through  the  data,  while  data  for 
emulsified  fuels  are  closed  symbols  with  dashed  or  dotted  lines.  Lines 
through  comparable  data  for  the  emulsified  and  neat  fuels  at  the  same 
condition  are  called  out  in  the  figures.  Drop  size  data  are  shown  as 
close  as  6.35  mm  (1/4  inch)  from  the  nozzle,  In  cases  where  the  spray 
density  was  not  too  great,  to  a  distance  of  31.35  mm  (1-1/4  inches)  or 
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56.35  mm  (2-1/4  Inches)  from  the  nozzle.  Data  are  not  shown  beyond  a 
point  where  excessive  noise  due  to  beam  steering  of  the  laser  beam 
prevented  reasonable  results. 


Atomization/Evaporation  Results  for  Emulsified  and  Neat  Jet  A 


The  effects  of  air  temperature  and  pressure  on  the  atomization  and 
evaporation  of  Jet  A  emulsions  and  neat  fuels  were  examined  using  nominal 
3.79  liter/hr  (1.00  gph)  and  11.4  liter/hr  (3.00  gph)  nozzles.  Drop  size 
data  for  the  11.4  liter/hr  nozzle  spraying  into  air  at  310  kPa  (3.06  atm) 
with  a  velocity  of  14.6  m/s  over  a  temperature  range  from  308  K  to  644  K 
are  shown  in  Figure  7.  Similar  data  for  the  3.79  liter/hr  nozzle 
spraying  into  higher  pressure  air  (448  kPa  or  4.42  atm)  moving  at  8.5  m/s 
are  shown  in  Figure  8.  In  Figure  8,  the  SMD's  have  been  normalized  so 
that  the  first  point  (high  temperature  cases)  or  average  of  first  two 
points  (lower  temperatures  cases)  is  equal  to  unity  to  show  the  relative 
evaporation  more  clearly.  The  initial  average  drop  sizes  in  Figure  8 
were  smaller  than  in  Figure  7  and  equal  to  about  31pm.  (The  smaller 
drop  sizes  are  typical  for  a  smaller  capacity  nozzle.) 

These  results  point  to  a  substantial  difference  in  the  evaporation 
characteristics  of  emulsified  jet  fuel  when  compared  with  neat  jet  fuel. 
When  the  air  temperature  is  low  (308  K) ,  evaporation  is  relatively 
unimportant  and  microexplosions  are  not  possible;  the  two  fuels  show  a 
similar  trend  in  drop-size  distribution,  indicated  by  the  Sauter  Mean 
Diameter  (SMD) .  These  data  and  the  more  extensive  results  presented  in 
Ref.  3  show  that  the  room  temperature  atomization  characteristics  of  the 
two  fuels  are  quite  similar  with  the  emulsified  fuel  usually  producing 
slightly  larger  fuel  drops  due  to  the  higher  viscosity.  This  fact  along 
with  the  observation  that  the  SMD's  near  the  nozzle  tend  to  be 
approximately  the  same  for  the  two  fuels  and  diverge  in  moving  away  from 
the  nozzle  for  the  higher  temperature  cases  indicate  that  the  differences 
in  SMD  are  due  to  differences  in  evaporation  rather  than  initial 
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FIGURE  8. 


atomization.  At  elevated  temperatures  there  is  a  systematic  difference 
between  two,  with  the  emulsified  fuel  always  showing  a  much  lower  average 
drop  size.  This  is  consistent  with  the  microexplosion  hypothesis,  but  it 
is  necessary  to  diverge  a  bit  and  try  to  explain  the  general  behavior  of 
all  the  data  in  Figures  7  and  8  to  better  understand  what  other  phenomena 
could  cause  the  same  results. 

First,  consider  that  these  drop-size  data  were  obtained  with  a 
laser-diffraction  technique  which  measures  drop  characteristics  along  its 
9  mm  diameter  line-of-sight,  and  that  the  measurement  method  can  affect 
the  results.  The  drop-size  data  in  Ref.  3  for  an  almost  nonevaporating 
spray  in  room  temperature  air  show  that  the  average  drop  size  sampled  by 
the  instrument  drops  rapidly  after  leaving  the  nozzle,  then  levels  out, 
then  Increases  slightly.  The  initial  reduction  in  SMD  near  the  nozzle  is 
due  to  at  least  three  phenomena.  First,  in  a  hollow  cone  nozzle  the 
majority  of  the  spray  volume  is  in  the  outer  part  of  the  cone  composed  of 
relatively  larger  drops,  while  a  smaller  amount  of  spray  made  up  of 
smaller  drops  fills  in  the  inner  core.  Near  the  nozzle  the  laser  beam 
intercepts  more  of  the  outer  cone  increasing  the  average  drop  size. 
Secondly,  secondary  atomization  may  occur  after  the  initial  breakup  which 
further  reduces  the  large  drops  as  distance  from  the  nozzle  Increases. 
Third,  for  cases  where  the  air  velocity  is  lower  than  the  initial  fuel 
velocity  as  In  Figures  7  and  8  and  Ref.  3,  the  small  drops  rapidly 
decelerate  to  the  air  velocity  increasing  their  relative  concentration  in 
the  sample  volume  and  decreasing  the  SMD.  The  differences  in  the  drag 
effects  on  large  and  small  drops  were  shown  dramatically  in  the  results 
of  the  computer  model  presented  earlier.  These  three  factors  contribute 
to  the  decrease  in  SMD  when  moving  away  from  the  nozzle  for  the  room 
temperature  case  where  evaporation  is  negligible. 

Now  consider  the  elevated  temperature  data  in  Figures  7  and  8  which 
Indicate  an  Increase  in  SMD  as  the  temperature  is  raised  and  the 
evaporation  rate  increased.  This  result  is  due  to  the  fact  that  although 
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all  drops  decrease  In  size  with  evaporation,  the  smaller  drops  evaporate 
much  more  quickly  and  the  resulting  average  drop  size  actually  increases. 
Thus  at  higher  air  temperatures  the  SMD  Increases  in  moving  away  from  the 
nozzle  relative  to  the  room  temperature  case.  This  is  demonstrated 
qualitatively  by  results  from  the  computer  model.  Data  in  Figures  7  and 
8  show  only  initial  trends  in  the  change  in  SMD  relative  to  the  total 
spray  lifetime.  The  computer  model  predicts  and  some  of  the  data 
indicate  that  the  SMD  initially  increases  and  then  levels  out  further 
into  the  evaporation  process.  Also  the  largest  SMD's  shown  are,  in  some 
cases,  biased  on  the  high  side  due  to  uncorrectable  beam  steering 
problems.  This  study  is  directed  at  the  relative  evaporation 
characteristics  of  emulsified  versus  neat  fuels  and  the  largest  SMD's 
reported  may  be  erroneously  high. 

Having  qualitatively  examined  the  phenomena  responsible  for  the 
trends  shown  for  neat  Jet-A  in  Figures  7  and  8,  it  is  possible  to 
consider  the  reasons  for  the  differences  between  the  emulsified  and  neat 
jet  fuel.  Since  the  two  fuels  show  a  similar  SMD  near  the  nozzle  and 
then  the  emulsified  fuel  shows  a  smaller  SMD  away  from  the  nozzle,  it  is 
entertaining  to  speculate  that  microexplosions  are  occurring  after  the 
emulsified  drops  heat  up  and  cause  the  larger  drops  to  subdivide  into 
smaller  drops.  However,  the  heat  of  vaporization  of  water  is  about  eight 
items  that  of  Jet  A  and  the  heat  capacity  of  water  is  also  higher  so  that 
the  water  may  be  absorbing  much  of  the  thermal  energy  available  and 
reducing  the  evaporation  rate  of  the  emulsified  fuel  so  that  the  rate  of 
increase  in  drop  size  is  reduced.  This  effect  must  be  balanced  against 
the  lower  boiling  point  of  water  which  would  tend  to  Increase  its 
evaporation  rate.  An  analysis  of  this  effect  will  be  performed  using  the 
computer  model  modified  to  handle  two  fuel  components. 

Further  tests  were  performed  with  emulsified  and  neat  Jet  A  at 
constant  air  temperature  with  varying  air  pressure.  Law  (Ref.  7)  has 
suggested  that  an  increase  in  air  pressure  at  constant  temperature  should 


favor  the  occurrence  of  microexplosions  because  the  steady  state  drop 
temperature  increases  with  air  pressure  as  was  shown  in  Figures  1  to  3. 
The  limit  of  superheat  of  water,  also  shown  in  Figures  1  to  3,  is  one  of 
the  criteria  used  to  predict  the  condition  for  the  onset  of 
microexplosions.  Spontaneous  nucleatlon  occurs  below  or  at  the  limit  of 
superheat  (about  88  to  90  percent  of  the  critical  temperature)  above 
which  fluids  will  always  film  boil.  Many  liquids  may  be  heated  well 
above  their  "boiling  points"  and  close  to  their  predicted  limit  of 
superheat  where  they  boil  explosively.  The  details  of  the  interface 
between  the  fluid  and  its  surroundings  usually  determine  where  the  fluid 
begins  to  boil  in  the  range  between  the  boiling  point  and  the  limit  of 
superheat.  Several  investigators  have  chosen  different  temperatures  to 
represent  the  condition  for  vaporization  of  the  water  and  drop 
disruption;  Birchley  and  Riley  (Ref.  8)  have  chosen  the  normal  boiling 
point,  Jacques  (Ref.  9)  and  Gollahalli  (Ref.  10)  use  that  temperature 
where  surface  tension  forces  are  exceeded  by  the  forces  due  to  the 
Internal  pressure  of  the  drop,  and  Law  (Ref.  7)  assumes  drop  disruption 
occurs  at  the  limit  of  superheat.  Measured  values  for  spontaneous 
nucleatlon  of  water  have  been  reported  by  Avedlslan  (Ref.  11). 

Thus,  although  there  is  disagreement  about  the  drop  temperature 
where  microexplosions  occur,  a  higher  drop  temperature,  as  allowed  by  an 
increase  in  air  pressure,  should  favor  a  divergence  in  drop  sizes  for 
emulsified  and  neat  fuels  if  drop  disruptions  are  occurring.  In  order  to 
examine  this  possibility,  tests  were  performed  at  589  K  at  pressures  of 
207  kPa  (2.04  atm)  and  310  kPa  (3.06  atm)  with  the  11.4  liter/hr  nozzle, 
and  at  547  K  with  pressures  from  165  kPa  (1.63  atm)  to  448  kPa  (4.42  atm) 
for  the  3.79  liter/hr  nozzle,  as  shown  in  Figures  9  and  10,  respectively. 
These  results  suggest  that  the  emulsified  and  neat  fuels  evaporate  in  a 
similar  manner  at  the  lower  pressures  but  tend  to  diverge  as  the  air 
pressure  is  increased  with  the  emulsified  fuels  producing  smaller  drops 
than  the  neat  fuels,  consistent  with  the  microexplosion  hypothesis. 
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FIGURE  10 


As  an  aside  to  the  comparison  of  neat  and  emulsified  fuel 
evaporation,  consider  the  effect  of  pressure  on  the  evaporation  of  neat 
fuel.  Figures  4  through  6  indicate  that  the  evaporation  constant  at 
quiescent  conditions  is  approximately  Independent  of  air  pressure  at 
temperatures  used  in  these  experiments,  while  Figures  9  and  10  show  that 
neat  fuels  evaporate  faster  as  the  pressure  is  increased  (assuming  that 
evaporation  rates  can  be  Inferred  from  the  rate  of  initial  increase  in 
SMD).  This  contradiction  can  be  partly  explained  by  three  factors. 
First,  as  the  air  pressure  Increases  the  initial  SMD  decreases  by  a 
factor  proportional  to  about  P"0*4,  and  the  smaller  drops  evaporate  much 
faster  according  to  the  d^-law.  Secondly,  the  computations  in  Figures  4 
through  6  are  for  the  quiescent  evaporation  constant  while  Eq.  3  shows 
that  the  convective  evaporation  constant  depends  on  pressure  through  the 
Reynolds  number.  Third,  the  transit  times  for  the  drops  to  a  fixed 
measurement  location  Increase  slightly  with  pressure  due  to  Increased 
drag  at  higher  pressures,  allowing  more  time  for  evaporation. 

Atomlzatlon/Evaporatlon  Results  for  Emulsified  and  Neat  Hexadecane 

In  order  to  make  these  arguments  more  quantitative,  it  was  necessary 
to  repeat  the  previous  experiments  with  a  single  component  fuel  whose 
evaporation  process  could  be  modelled  with  greater  certainty.  Three 
single  component  fuels  were  obtained:  n-hexadecane,  n-dodecane,  and  n- 
octane,  but  experiments  have  been  performed  only  with  the  n-hexadecane  at 
this  time.  The  evaporation  of  hexadecane  drops  was  investigated  at  fixed 
air  pressure  and  varying  air  temperature  and  visa  versa.  At  a  fixed 
pressure  of  4.42  atm  (448  kPa)  the  Sauter  mean  diameter  (SMD)  of  the 
spray  varied  with  distance  from  the  nozzle  as  shown  in  Figure  11.  A 
comparison  of  the  size  distribution  of  emulsified  and  neat  fuels  at  506  K 
and  547  K  indicates  no  difference  within  the  precision  of  the 
measurement.  (Data  for  the  emulsion  at  506  K  were  essentially  identical 
with  those  for  the  emulsion  at  547  X  and  were  omitted  in  Figure  11.)  At 
589  K  the  emulsified  fuel  produces  smaller  drops  and  the  difference 
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FIGURE  11.  COMPARISON  OF  RELATIVE  SMD  FOR  EMULSIFIED  AND  NEAT  HEXADECANE 
AT  CONSTANT  PRESSURE  (4.42  atm,  448  kPa)  AND  VARIOUS  TEMPERATURES 
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increases  at  the  higher  temperatures.  The  corresponding  predicted 
maximum  steady-state  surface  temperatures  of  the  drops  as  computed  by  the 
model  are  shown  in  Table  3.  At  a  fixed  temperature  of  644  K  the  drop 
sizes  of  neat  and  emulsified  fuels  were  studied  as  a  function  of  air 
pressure  with  the  results  as  shown  in  Figures  12-14.  Examination  of 
Figure  14  indicates  that  at  the  lower  pressures  of  165  kPa  (1.63  atm), 
241  kPa  (2.38  atm),  and  345  kPa  (3.40  atm),  the  drop  sizes  in  sprays  of 
neat  and  emulsified  fuels  are  about  the  same,  but  that  at  higher 
pressures  the  sprays  of  emulsified  fuels  have  smaller  SMD's  during  the 
heat  up  and  evaporation  process.  Table  4  presents  the  computed  droplet 
temperatures  for  these  experiments. 

Table  3.  Computer  Predicted  Steady  State  Drop  Surface  Temperatures 
for  Hexadecane  at  Conditions  Corresponding  to  Figure  11 
<pair  "  448  kPa  (4.42  atm)) 

Air  Temperature  Drop  Temperature 
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00 

506 

488 

547 

503 

589 

520 

644 

535 

700 
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Table  4.  Computer  Predicted  Steady  State  Drop  Surface  Temperature 
for  Hexadecane  at  Conditions  Corresponding  to  Figures  12-14 

(T.ir  -  M4  K) 
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Air  Pressure 
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Drop  Temperature 
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165 

1.63 

491 

241 

2.38 

502 

345 

3.40 

512 

448 

4.42 

520 

586 

5.78 

529 
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FIGURE  14.  RESULTS  OF  FIGURES  12  AND  13  COMPARED 


These  treads  correspond  qualitatively  to  those  expected  if  the 
emulsified  fuel  were  undergoing  microexplosions  when  the  temperature  and/ 
or  pressure  were  raised  to  a  point  where  the  drop  temperature  exceeded 
the  spontaneous  nucleation  point  of  water.  Below  the  spontaneous 
nucleation  point  the  spray  drop  sizes  should  be  similar  for  the  two 
fuels,  hut  as  the  drop  temperature  exceeds  that  condition  the  emulsified 
drop  may  be  shattered  into  many  smaller  drops  producing  a  smaller  SMD 
than  the  neat  fuel  as  distances  from  the  nozzle  increase. 

In  order  to  make  this  argument  quantitative,  it  is  necessary  to  know 
the  drop  temperature  as  a  function  of  the  air  temperature  and  pressure. 
Presently  there  are  no  techniques  to  measure  drop  temperatures  in  a 
polydisperse  spray  where  the  small  drops  heat  up  rapidly  while  the  large 
ones  are  much  slower.  This  makes  it  necessary  to  predict  the  drop 
temperatures  using  the  computer  model  to  balance  the  heat  conducted  to 
the  drop  from  the  air  with  the  heat  used  to  vaporize  the  fuel.  Some 
results  from  the  model  described  earlier  are  shown  in  Tables  3  and  4 
which  indicate  the  Increase  in  steady-state  drop  temperature  with 
increasing  air  temperature  or  pressure.  By  comparing  those  tables  with 
Figures  11  and  14,  it  may  be  seen  that  the  drop  size  distributions  for 
the  emulsified  and  neat  fuels  diverge  at  a  predicted  steady  state  drop 
temperature  of  about  520  K.  This  may  be  compared  with  the  limit  of 
superheat  of  water  of  about  540  K  (Ref.  12).  The  divergence  of  drop  size 
distributions  for  emulsified  and  neat  fuels  at  a  drop  temperature 
somewhat  below  the  limit  of  superheat  is  consistent  with  the 
microexplosion  hypothesis  (Ref.  11).  As  mentioned  previously,  the  limit 
of  superheat  is  an  upper  limit  where  drop  fragmentation  must  occur  if 
water  is  present,  while  the  probable  temperature  for  spontaneous 
nucleation  of  water,  and  hence,  microexplosions  to  occur  is  bounded  by 
the  normal  boiling  point  (corrected  for  pressure)  on  the  low  end  (~373 
K)  and  the  limit  of  superheat  ( ~ 540  K)  on  the  upper  end.  However, 
results  from  the  model,  which  also  predicts  heat  up  times,  would  indicate 


only  the  smaller  diameter  drops  have  reached  the  steady-state  temperature 
at  positions  where  the  drop  size  distributions  diverge. 

Table  5  represents  the  computer  predicted  drop  temperatures,  (FTEMP) 
and  the  other  drop  Information  as  explained  previously  for  Table  2,  for  a 
location  (X  DROP)  26.4  mm  from  the  nozzle  for  air  conditions  of  700  K  and 
448  kPa  (4.42  atm).  At  these  conditions,  the  results  from  Figure  11 
Indicate  a  considerable  divergence  between  the  emulsified  and  neat  fuel 
SMD's.  Although  the  steady-state  temperature  is  close  to  the  superheat 
limit  of  water  (540  K)  at  these  conditions,  only  the  drops  of  Initial 
size  smaller  than  28.0  pm  exceed  (approximately)  500  K  and  they  represent 
about  54  percent  of  the  original  liquid  volume  of  the  spray.  The  drops 
equal  to  or  larger  than  177.6  pm  represent  less  than  three  percent  of  the 
liquid  volume  of  the  Initial  spray  and  can  be  Ignored.  Thus,  although 
the  steady-state  temperature  Is  about  535  K,  most  of  the  drops  are 
predicted  to  reach  that  temperature  only  late  In  their  lifetime,  and 
larger  drops  have  only  begun  to  heat  up;  e.g.,  the  14lp  m  drop  has 
Increased  In  "uniformly  mixed"  temperature  from  311  K  to  324K.  Peak 
surface  temperatures  in  the  larger  drops  are  higher,  of  course,  than  the 
uniformly  .nixed  temperature  would  imply,  but  the  result,  nevertheless,  Is 
that  the  larger  drops  in  the  spray  are,  for  the  conditions  of  these 
experiments,  much  lower  in  temperature  than  the  steady-state  temperature. 
The  variation  of  drop  temperature  with  size  coupled  with  the  uncertainty 
of  the  temperature  where  spontaneous  nucleatlon  of  the  water  microdrops 
should  occur  (as  discussed  earlier  in  the  report)  makes  an  exact 
condition  where  microexplosions  should  occur  a  meaningless  exercise. 
What  the  model  results  do  reveal  Is  that  the  steady  state  drop 
temperature  does  increase  with  Increasing  air  temperature  and  pressure  to 
a  range  where  spontaneous  nucleatlon  probably  would  occur  If  the  water 
has  not  completely  evaporated  before  reaching  that  condition.  The  model 
also  indicates  a  considerable  amount  of  the  liquid  spray  volume  which  Is 
contained  In  the  larger  drops  is  not  at  a  temperature  where 
microexplosions  are  expected. 


Table  5.  Ccaputer  Predicted  Spray  Characteristics  for  a  Condition  of 
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The  computer  model  Is  instructive  in  predicting  and  understanding 
some  of  the  phenomena  observed  experimentally,  but  it  is  not  sufficiently 
complete  to  compare  precisely  with  the  experimental  results.  The 
limitations  of  the  current  version  of  the  computer  code  were  mentioned 
earlier  in  this  report  and  include:  (1)  no  provision  for  inhomogenlties 
in  drop  size  distribution  along  the  line-of-sight ,  (2)  secondary 
atomization  of  drops  near  the  nozzle  is  ignored,  and  (3)  the  weighting 
factor  associated  with  the  laser  beam  location  and  size  has  been  omitted. 
Because  of  the  first  two  omissions,  the  model  predicts  an  increase  in  SMD 
moving  away  from  the  nozzle  while  the  experimental  results  show  an 
initial  decrease  except  at  the  highest  temperatures.  An  effort  was  made 
to  compare  the  trends  for  experimental  and  predicted  SMD's  by  adjusting 
the  initial  Rosin-Rammler  parameters  (x  and  N)  so  that  the  predicted  x 

and  N  at  11.4  mm  and  16.4  mm  were  similar  to  the  measured  values  with  the 

/ 

results  as  shown  in  Figure  15.  Experimentally  measured  SMD's  at  6.4  mm 
were  usually  larger  than  at  11.4  mm  while  the  predicted  SMD's  were 
significantly  smaller.  The  fact  that  the  atomization  process  was  not 
complete  at  the  6.4  mm  location  as  well  as  the  other  omissions  from  the 
model  just  mentioned  probably  account  for  that  initial  discrepancy.  It 
should  be  noted  that  the  upstream  edge  of  the  laser  beam  is  only  2  mm 
from  the  nozzle  at  the  6.4  mm  sampling  location.  There  is  a  considerable 
discrepancy  at  the  21.4  mm  location  for  the  700  K  case.  The  discrepancy 
is  probably  partly  due  to  an  overestimate  of  the  measured  SMD  due  to  beam 
steering  problems,  and  points  out  the  need  to  measure  the  amount  of  fuel 
vaporized  at  different  locations  to  further  check  out  the  model.  A  more 
exact  comparison,  of  experimental  and  model  results  will  require  some 
additional  sections  to  be  added  to  the  model  to  correct  the  deficiencies 
cited. 

Drop  Size  Measurements  in  Flames 

Drop  size  measurements  in  burning  sprays  have  been  attempted  in 
several  configurations  including  cases  where  the  flame  was  stabilized  in 


a  region  downstream  from  the  nozzle.  Experiments  to  date  have  resulted 
in  inadequate  signal-to-noise  ratios  to  determine  reasonable  drop  size 
distributions.  Generally,  two  types  of  problems  are  present  when  making 
size  measurements  with  laser-diffraction  instrumentation  in  burning 
sprays.  First,  the  flame  produces  high  levels  of  visible  and  near- 
infrared  radiation  (in  the  response  range  of  the  detectors)  which  is 
orders  of  magnitude  higher  in  Intensity  than  the  desired  scattered 
radiation  from  the  HeNe  laser.  Secondly,  the  severe  gradients  in 
temperature  and  air/fuel  vapor  concentrations  result  in  serious  beam 
steering  problems  causing  the  laser  beam  to  be  directed  onto  areas  of  the 
detector  other  than  the  central  ring.  This  produces  signals  which  appear 
to  be  due  to  large  drops  scattering  light  onto  the  inner  detector  rings 
and  thus  represents  a  false  signal. 

The  three  modifications  described  in  the  EXPERIMENTAL  APPARATUS 
section  were  sufficient  to  discriminate  against  the  flame  radiation  for 
the  cases  examined.  However,  the  severe  refractive  index  gradients 
produced  so  much  noise  on  the  inner  detector  rings  that  the  drop  size 
data  which  were  computed  were  judged  to  be  unreasonable.  Further 
consideration  is  being  given  to  ways  to  circumvent  this  problem. 


CONCLUSIONS 


Techniques  have  been  developed  which  allow  for  the  determination  of 
drop-size  data  in  evaporating  fuel  sprays  over  a  range  of  elevated 
pressures  and  temperatures  important  to  gas-turbine  combustion. 

It  has  been  determined  that  there  are  substantial  differences  in  the 
evaporation  characteristics  of  emulsified  fuels  when  compared  with  neat 
fuels.  The  spray  droplets  of  the  evaporating  emulsified  fuel  are,  on  the 
average,  significantly  smaller  than  the  droplets  of  neat  fuel  at 
corresponding  elevated  temperatur e/pres sure  conditions.  The  differences 
are  enhanced  as  the  pressure  increases  for  the  range  examined  (to  about 
5.8  atm  (590  kPa)).  These  findings  are  consistent  with  the 
microexplosion  hypothesis  but  a  thorough  analysis  of  other  possible 
effects  has  not  been  completed. 

Limited  experimental  data  would  suggest  that  the  size  of  drops  at 
initial  atomization  depends  strongly  on  the  air  pressure,  with  SMD 


RECOMMENDATIONS  FOR  NEXT  TEAR 


Results  already  described  have  Indicated  significant  differences  in 
drop  sizes  for  evaporating  sprays  of  emulsified  and  neat  fuels.  In  order 
to  understand  the  reasons  for  these  differences,  the  following  work  Is 
proposed  for  the  next  year  of  this  program,  although  there  may  be 
Insufficient  time  to  pursue  all  of  these  suggestions. 

1.  Conduct  similar  atomization/evaporation  experiments  with  single 
component  fuels  of  lower  boiling  point  than  n-hexadecane.  Fuels 
such  as  n-octane/water  emulsions  should  not  exhibit  microexplosions 
but  would  still  be  subject  to  the  "heat-sink"  effect  of  water  on  the 
evaporation  process.  This  will  help  to  separate  the  microexplosion 
effect  on  drop  size  from  the  reduced  evaporation  rate  due  to  the 
high  heat  of  vaporization  of  water. 

2.  Develop  method(s)  for  measuring  fuel  vapor  concentration  In  sprays. 
Both  optical  and  probe  techniques  have  been  used  by  other  workers  to 
measure  fuel/alr  ratios  in  sprays  and  to  differentiate  between  the 
liquid  and  vapor  phases  of  the  fuel.  These  would  be  useful  in 
studying  the  evaporation  of  emulsified  fuels,  as  microexplosions 
would  enhance  the  evaporation  rates  while  the  "heat-sink"  effect  of 
water  would  slow  down  evaporation  rates. 

3.  Modify  computer  model  for  2-component  fuels  with  various  assumptions 
about  relative  evaporation  rates.  The  effect  of  water 
emulsification  with  a  single  component  fuel  could  be  studied  from  a 
theoretical  view  to  compare  with  the  experimental  results. 

4.  Modify  computer  model  tc  Include  the  inhomogeniety  of  spray,  and 
laser  beam  location  and  intensity  profile.  This  would  lead  to  a 
better  capability  to  compare  experimental  and  computer  predicted 
results. 


Modify  computer  model  to  predict  the  concentration  of  fuel  vapor  as 
a  function  of  axial  distance  from  the  nozzle.  If  the  measurement 
techniques  described  in  item  (2.)  are  developed,  they  could  be 
compared  with  the  predicted  results. 
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APPENDIX  A 


FUEL  SPRAY  COMPUTER  MODEL  -  MATHEMATICS 


For  purposes  of  discussion,  the  model  Is  divided  into  three  parts 
called  (1)  Thermodynamics,  (2)  Aerodynamics,  and  (3)  Drop-Size 
Distribution.  Each  section  Is  discussed  in  order.  An  overall  flow  chart 
Is  shown  In  Figure  A-l. 


The  mods 


This  part  of  the  model  computes  the  heat-up  of  the  drop,  the  final 
steady-state  temperature,  and  the  fuel  and  air  properties  necessary  for 
those  calculations.  This  part  of  the  model  is  identical  to  that 
described  by  Chin  and  Lefebvre  (Ref.  A1  and  A2)  except  that  some  of  the 
fuel  property  data  were  taken  from  the  American  Petroleum  Institute  Data 
Book  (Ref.  A3).  The  results  from  this  model  are  very  similar  to  those 
presented  by  Chin  and  Lefebvre  when  the  model  is  used  for  quiescent  fuel 
sprays.  The  integration  of  this  part  of  the  model  with  the  other 
sections  to  be  described  allows  the  examination  of  sprays  in  situations 
where  both  spray  and  air  are  moving  with  a  nonzero  relative  velocity. 
The  mathematics  are  outlined  below,  but  a  more  detailed  derivation  is 
given  by  Chin  and  Lefebvre  (Ref.  A1  and  A2)  which  in  turn  is  based  on  the 
theories  described  by  Spalding  (Ref.  A4). 

The  steady-s tate  temperature  Is  determined  as  follows.  The  mass 
transfer  rate,  m,  from  a  drop  is  gi/en  by, 


mp  =  2ttD  (^•g/cp>g)  In  (1+B) 


where  D  is  the  diameter,  kg  is  the  thermal  conductivity  of  the  gas,  Cp,g 
is  the  specific  heat  at  constant  pressure,  and  B  is  either  B^,  the  mass 
diffusion  transfer  number  if  mass  diffusion  is  controlling,  or  Bj,  the 


FIGURE  A-l  FLOW  CHART  OF  SPRAY  MODEL 


thermal  diffusion  transfer  number  If  heat  transfer  is  controlling  and  the 
Lewis  number  is  assumed  to  be  unity*  is  defined  as 

BM  *  (Yf,s  -  YF>»)  /  (1  -  Yf,s)  (A2) 

where  Yp>s  and  YF>  are  the  fuel  mass  fractions  at  the  drop  surface  and 
in  the  ambient  air,  respectively,  and  Yp(S  is 

YF,s  =  Pf,s  Mf  /  (PF,s  MF  +  (P  -  PF>S)  Ma)  (A3) 

where  Pf,s  Is  the  fuel  vapor  pressure  at  the  drop  surface,  P  is  the 
ambient  pressure,  and  Mp  and  Ma  are  the  molecular  weights  of  fuel  and 
air,  respectively.  The  vapor  pressure  was  calculated  using  the  method 
recommended  by  the  American  Petroleum  Institute  Data  Book  (Ref.  A3) 
called  Procedure  5A1.10,  "Vapor  Pressures  of  Pure  Hydrocarbons." 

The  thermal  diffusion  transfer  number  Bx  for  an  evaporating  drop  is 

bT  -  cp>g  (T„  -  Ts)  /  L  (A4) 

where  Too  and  Ta  are  the  ambient  air  temperature  and  drop  surface 
temperature  and  L  is  the  latent  heat  of  fuel  vaporization  corrected  from 
the  normal  boiling  point  to  the  actual  surface  temperature. 

The  accuracy  of  Equation  (A_l)  is  very  dependent  on  the  choice  of 
values  for  kg  and  Cp^.  As  recommended  by  Chin  and  Lefebvre  (Ref.  Al)  a 
reference  temperature  was  chosen  as  the  drop  surface  temperature  plus  1/3 
of  the  difference  between  ambient  air  and  surface  temperatures. 
Similarly  a  reference  value  of  the  fuel  vapor  mass  fraction  was  taken  as 
the  value  just  outside  the  surface  plus  1/3  of  the  difference  between  the 
value  at  infinity  and  just  outside  the  surface.  Using  these  reference 
conditions  denoted  by  "r"  the  specific  heat  of  the  gas  is  given  by, 


Cp.g  -  YA,r  (Cp,A  at  Tr)  +  Yp  >r  (cp>v  at  Tr) 


(AS) 


and  the  thermal  conductivity  by, 

kg  =  YA,r  (^A  at  Tr)  +  Yp>r  (ky  at  Tr)  (A6) 


At  steady-state  conditions,  Bh  -  B-j  or  from  Equations  (A2)  and  (A4) 


YF,s  -  YF,®  _  cp,g  (T»  -  Ts>st) 
l-YF,s  L 


(A7) 


where  Ts>st  is  the  desired  surface  temperature  under  steady-state 
conditions.  Assuming  that  the  fuel  vapor  concentration  in  the  ambient 
air  is  zero  (Yp>00  -  0)  and  substituting  for  Yp|g  from  Equation  (A3)  then 
(A7)  becomes 

P  Mp  L 

-  _  —  -  _  1  =  n  (A8) 

PF,s  MA  cpjg  (T®  -  Ts,st) 


In  solving  Equation  (A8)  to  determine  the  steady-state  surface 
temperature  Ts^stf  some  of  the  variables  can  be  specified,  P,  Mp,  Ma,  and 
Too  but  the  remaining  three  variables  are  functions  of  temperature,  with 
L  being  related  to  the  heat  of  vaporization  at  the  normal  boiling  point 


t'T.bn  by. 


L  *  LT,bn 


[1  0.38 

(Tcr  “  Ys)  /  (Tcr  "  Tbn)  J 


(A9) 


where  Tcr  is  the  critical  temperature  and  Tt,n  is  the  normal  boiling  point 
of  the  fuel.  In  order  to  solve  (A8)  it  is  necessary  to  first  assume  a 
value  for  T8>st.  That  temperature  is  used  to  calculate  a  reference 
temperature  Tr  which  in  turn  is  used  to  calculate  those  quantities 
specified  in  Equations  (A5)  and  (A6)  used  to  calculate  cp>g  and  kg,  and 
to  calculate  L  as  specified  in  (^9).  If  the  assumed  value  of  Ts>st  is 
too  low,  the  left  side  of  Equation  (A8)  will  be  positive.  An  iterative 
procedure  must  be  used  to  continue  specifying  TSjSt  until  Equation  (A8) 
is  satisfied  within  the  desired  accuracy.  That  procedure  results  in  a 


1 

1 

determination  of  the  steady  state  surface  temperature  Tg>gt  and  the 
transfer  number  B*Bm**Bt.  The  steady  state  quiescent  evaporation  constant 
Xq>gt  can  then  be  determined  from, 

Aq,st  =  8  kg  In  (1+B)  /  Pp  Cp,g  (A1°) 

where  is  the  fuel  density.  Frossling  (Ref.  A5)  has  shown  that 
convective  effects  can  be  accounted  for  in  the  case  where  heat  transfer 
rates  are  controlling  by, 

xc,st  -  xq,st  C1  ♦  °-276  ReD°*5  Prg0-53)  (All) 

where  Xc>gt  is  the  steady-state  evaporation  constant  corrected  for 
convective  effects,  Prg  is  the  Prandtl  number  for  the  gas,  and  Rep  is  the 
Reynolds  number  using  the  relative  velocity  between  the  drop  and  the  gas. 

This  velocity  should  include  a  fluctuating  component  in  the  case  of 
turbulent  flow.  Turbulence  intensities  were  not  measured  in  these 
experiments,  but  a  value  of  20  percent  of  the  axial  velocity  was  used  as 
an  estimate.  The  effective  velocity  used  to  calculate  Rep  for  the 
evaporation  calculation  was  taken  as  the  sum  of  the  fluctuating  component 
(20  percent  of  the  air  velocity)  plus  the  difference  between  the  drop  and 
air  velocities.  These  computations  specify  the  steady-state  properties 
of  the  drop. 

The  evaporation  which  occurs  during  the  drop  heat-up  period  is 
significant  in  many  practical  situations,  with  the  drop  being  completely 
evaporated  before  reaching  steady-state  conditions  is  some  instances. 

Chin  an<*  Lefebvre  (Ref.  A2)  have  predicted  that  high  air  pressures  and 
convective  effects,  which  are  often  absent  in  laboratory  experiments, 
both  tend  to  increase  the  relative  Importance  of  the  heat-up  period 
relative  to  the  steady-state  phase  of  drop  evaporation.  Because  the 
experiments  conducted  in  this  study  were  of  real  sprays  at  high  pressures 
and  relative  velocities,  it  was  necessary  to  model  the  heat-up  period  in 
detail. 


A-5 


Chin  and  Lefebvre  (Ref.  A2)  have  shown  that  the  rate  of  change  of 
drop  surface  temperature  Is  given  by, 

dTs  /  dt  -  (iF  L  /  Cp,F  m)  {  (Bt  /  Bm)  -  1  )  (A12) 

where  mp  is  specified  by  Equation  (Al)  with  B*Bm  during  the  heat-up 
period,  and  m  is  the  drop  mass, 

m  =  {w/ 6)  PF  D3  (A13) 

The  change  of  drop  size  with  time  is  given  by  (Ref.  A2), 

dD  /  dt  =  -4  kg  In  (1+BM)  /  *>p  cp>g  D  (A14) 

The  drop  temperature  asymptotically  approaches  the  steady-state 
temperature,  and  it  was  assumed  that  when  the  temperature  had  risen  to  95 
percent  of  the  difference  between  the  initial  fuel  temperature  and  the 
steady-state  temperature,  that  the  heat-up  calculations  could  be 
terminated  and  steady-state  properties  used. 

Aerodynamics 

The  approach  used  was  a  much  more  simplified  version  of  the  one 
developed  at  the  University  of  Sheffield  and  described  in  several  papers, 
e.g.,  Boyson  and  Swithenbank  (Ref.  A6).  A  cylindrical  coordinate  system 
was  used  where  x  is  the  axial  distance,  y  the  radial  distance,  and  z  the 
angular  position,  with  corresponding  velocities  in  the  axial  (u),  radial 
(v),  and  tangential  (w)  directions.  For  this  work,  the  tangential 
velocities  were  assumed  zero. 


The  equations  of  motion  of  a  particle  (drop)  neglecting  ell  forces 
except  drag,  are, 


F  (vp  -  vj 


(A16) 


yp 


Vp  Wp 

W  =  - -  -  F  (Wp  -  *„)  (A17) 

yP 


where  the  "p"  subscripts  refer  to  the  particle  and  the  ”<»**  subscript 
refers  to  the  free  stream,  and  F  is  given  by, 


F  = 


(18  yg  /  *>p 


2 

Dp 


)  (CD  Re  /  24) 


(A18) 


where  is  the  gas  viscosity,  Pp  is  the  particle  density,  Dp  is  the 
particle  diameter,  Cd  is  the  drag  coefficient,  and  Re  is  the  Reynolds 
number  defined  as. 


Re 


DP  Pg 


(A19) 


The  drag  coefficient,  Ct),  is  given  by  (Ref.  A7), 


CD  =  27  Re'0’84 

0  Re  - 

80 

(A20) 

CD  =  0.271  Re0-217 

80  <  Re  * 

104 

(A21) 

CD  =  2 

104  >  Re 

(A22) 

The  equations  of  motion  (A15  to  A17)  are  solved  numerically  using  a  step 
size  of  1  to  10  ps,  and  the  equations  of  trajectory  are  solved  in  a 
similar  manner, 
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These  equations  describe  the  trajectory  of  a  particle  in  a  gas  stream. 
For  each  iteration  in  the  trajectory  calculation,  the  drop  size  and 
temperature  are  updated  using  the  procedure  described  in  the 
Thermodynamics  section. 

At  the  end  of  each  iteration  through  the  trajectory  calculations, 
the  current  calculated  axial  position  of  the  particle  is  checked  against 

a  target  value.  When  the  particle  reaches  that  target  value  for  axial 

/ 

position,  the  position,  transit  time,  size,  temperature,  and  other 
parameters  for  that  drop  are  frozen  and  the  calculations  are  repeated  for 
the  next  larger  drop  size.  After  all  drops  reach  a  targeted  axial 
position,  corresponding  to  a  measurement  location,  the  drop  size 
distribution  is  calculated,  Important  data  are  printed  out,  the  target 
position  is  moved  downstream,  and  the  calculations  are  repeated. 


Prop  Size  Distribution 

The  computer  model  makes  all  calculations  on  individual  drops  of 
certain  sizes,  but  the  instrumentation  used  to  measure  the  spray  works 
best  in  characterizing  the  drop  sizes  by  two  parameters  which  specify  a 
size  and  width  of  the  distribution.  Three  distributions  may  be  selected: 
normal,  log-normal,  and  Rosin-Rammler.  The  Rosln-Rammler  distribution 
has  provided  the  best  fit  of  the  sprays  studied  and  has  been  used 
throughout  this  program.  In  order  to  compare  the  predictions  from  the 
computer  model  with  the  experimental  results,  it  is  necessary  to  convert 
the  computer  predicted  drop  sizes  at  a  given  location  into  an  equivalent 
set  of  two  parameters  specifying  the  Rosln-Rammler  distribution. 


If  R  represents  the  fraction  of  the  liquid  being  sprayed  contained 
in  drops  larger  than  diameter  D,  then  for  the  Rosin-Rammler  distribution 
(Ref.  A8), 

R  *  exp  (- CD  /  x)N  )  (A26) 


where  x  represents  a  size,  and  N  specifies  the  width  of  the  distribution. 
It  is  convenient  to  specify  a  spray  by  a  single  parameter  representing  an 
"average”  drop  size.  However,  a  straight  numerical  average  is  heavily 
weighted  towards  the  smallest  drops  which  are  extremely  plentiful  but 
which  contain  a  very  small  fraction  of  the  total  volume  of  the  liquid. 
For  combustion  processes,  the  surf  ace-  to-volume  ratio  la  important  to 
evaporation,  so  an  average  drop  size  may  be  chosen  which  has  a  surface- 
to-volume  ratio  representative  of  the  actual  spray,  and  such  an  average 
is  called  the  surface-volume  mean  diameter  or  Sauter  mean  diameter  (SMD). 
The  SMD  is  related  to  the  Rosin-Rammler  parameters  by  (Ref.  A8), 

SMD  =  x  /  T  (1-1/N)  N>1  (*27) 

Specifying  x  and  N  allows  the  calculation  of  the  SMD,  or  specifying  SMD 
and  N  determines  x. 

An  initial  set  of  drop  sizes  distributed  approximately  exponentially 
In  size  (i.e.,  evenly  spaced  when  plotted  as  In  D)  and  covering  the  range 
of  about  6ym  to  560ym  (corresponding  to  the  response  range  of  the  drop 
sizing  instrument)  with  16  to  21  different  sizes  is  used  in  the  model. 
Assuming  size  classes  bounded  by  these  drops  as  end  points,  the  initial 
fraction  of  liquid  contained  in  drops  larger  than  a  certain  size  class  is 
given  by  Equation  (A26). 


Ri  =  exp  (-(D0>i  /  x)N) 
R0  5  1 
Rn+1  =  0 
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ISiSn 


(A28) 


where  ”n"  is  the  number  of  drops,  D0  is  the  initial  drop  diameter,  and 
the  initial  fraction  of  liquid  in  each  size  class  is, 

Fo,i  *  Ri  “  Ri+1  o < i < n  (A29) 


After  evaporation  begins,  the  small  drops  begin  evaporating  quickly 
while  the  large  ones  evaporate  slowly,  changing  the  fraction  of  liquid  in 
the  different  size  classes  unevenly.  Denoting  the  smallest  nonevaporated 
drop  size  by  ”k”,  the  fraction  of  liquid  remaining  in  any  size  class  is, 

Fo,i  {(Di  /  D0>i)3  /  Up  i  ♦  (Di+1  /  D0ji+1)3  /  up>i+1  j/2 

=  -  - - - - - - -  (A30) 

F°»j  ^  ^o,j)^  /  Up.i  +  (Dj+i  /  D0>j+j)3  /  Up f j ♦  i  f/2 

j=k 

where 

D0  5  Dl 
UP.O  5  Up,l 
°n+l  =  Dn 
“p.n+l  5  up,n 

and  up>i  is  the  axial  velocity  of  the  ith  drop,  D*  is  the  Instantaneous 
drop  size,  and  D0|£  is  the  initial  drop  size  of  the  1th  drop.  Each  size 
class  is  characterized  by  the  drops  making  up  the  end  points,  and  the 
initial  fraction  of  liquid  in  a  size  class  is  modified  by  the  average 
loss  of  volume  of  the  end  point  drops.  The  velocity  term  is  added  to 
account  for  the  fact  that  as  the  drops  slow  down  they  Increase  their 
relative  concentration,  and  thus  their  weighting  factor,  in  the  sample 
volume. 


k-SiSn 

(oikin) 


The  cumulative  fraction  of  liquid  in  all  size  classes  including  and 
larger  than  the  ith  class  is, 


k<i<n 


(A31) 


This  instantaneous  value  of  R^  is  related  to  the  drop  sizes  and 
Rosln-Rammler  parameters  by  Equation  (A28)  with  the  initial  drop  size 
D0>1  replaced  by  the  Instantaneous  value  D^, 


Ri  =  exp  C-(Di  /  x)N)  k  <  i  <  n 


(A32) 


Taking  the  natural  logarithm  of  both  sides  twice  and  excluding  i*k 
(Ri«l),  Equation  (A32)  becomes, 

lne(lne(l  /  Ri))  =  N  lne  Di  -  N  lne  x  k+1  £  i  <  n  (A33) 

This  has  the  fora  of  the  equation  of  a  straight  line,  y-mx+b,  if  the 
following  definitions  are  used, 

y  -  lne(lne(l/R) ) 

m  -  N 

x  *  lne  D 
b  ■  -N  lne  x 

Thus,  by  determining  a  least  squares  fit  of  the  straight  line 
through  the  data  lne(lne(  1/Ri ) )  versus  lneDi,  the  Rosln-Rammler 
parameters  are  given  by, 

N  =  m  (the  slope)  (A35) 
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and 


x  =  exp  (-b  /  N) 


A  standard  routine  is  used  to  perform  the  least  squares  operation  and  the 
Ro8in-Rammler  parameters  and  the  SMD  (from  Equation  (A27))  are  determined 
at  each  target  value  of  the  axial  location  corresponding  to  the  position 
where  experimental  data  are  obtained.  Although  the  initial  distribution 
at  the  nozzle  is  an  ideal  Rosin-Rammler  distribution,  the  distribution 
downstream  does  not  correspond  exactly  to  the  Rosin-Rammler  equation  due 
to  the  different  evaporation  rates  for  the  different  sized  drops,  and  the 
degree  of  fit  is  determined  by  the  correlation  coefficient  of  the 
straight  line  through  the  computed  data. 
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APPENDIX  B 


FUEL  SPRAY  COMPUTER  MODEL  -  LISTING 


&DROPP  T= 00004  IS  ON  CR  TR  USING  00068  BLKS  R=000u 


000) 
0002 
00  03 
0004 
0005 
0006 
000? 
0008 
0009 
0010 
0011 
00  12 
0013 
0014 
0015 
0016 
0017 
0018 
0  019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0  033 
0  034 
0035 
0036 
0  037 
0038 
0039 
004  0 
0041 
0042 
0043 
0  044 
0045 
0046 
0047 
0043 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
0058 


FTN4 ,  L ,  Q 

PROGRAM  DROPP 

C  THIS  DROP  EVAPORATION  MODEL  INCORPORATES  MUCH  OF  THE  HEAT  TRANSFER 
C  SUGGESTED  BY  CHIN  AND  LEFEBVRE  IN  AIAA-S2-H76  -  8/3  0/82 

D I MENS  ION  D< 23>,DI AM  0<  23  > ,  T I ME<  23  > , XP<  23  > , YP< 23  » , ZP< 23  > , UPC  23  > , 
C  VP<  23  > ,  WP<  23  > , FF  0<  24  ) ,  FF<  24  >  ,  R<  24  > ,  FTEMP<  23  > ,  P.E-;  23  > ,  NU<  23  > , 

C  F< 24 > , CB< 23  > , CD<  23  > , XDATA< 23  > , YDATA< 23), ZDAT  A< 2  3  > , DM  IK  23  > , 

C  XK<  23  > , REM<  25  )  , 

&  TSD<  23  ) ,  LAMBHU<  23  ) ,  LAMB  A'.'  23  > 

REAL  L ,  LAMBST ,  NU ,  MUA ,  N ,  NNEW ,  M ,  LTBN ,  MWh  ,  MWF  ,  LEFT  ,  K.h  ,  K  V  ,  KG  ,  MUG  ,  MU 
&  LAMBHU, LAMBA, MASSD, MOOT 
C  DROP  DIAMETERS  <M> 

DATA  D I  AM  0/5 .8, 7. 2, 9. 0,1  1  . 4 , 14,4, 18,4,23.6,30.3,33.9,50. 1 , 

C  64.6,84.2,112.8,160.3,261 .6,563.9/ 

NDROP  =  16 
DO  10  1=1, NDROP 
)  0  DIAMCK  I  >  =  DIAM0<  I  )  /  1E6 
C  ALTERNATE  SET  OF  DROP  DIAMETERS 
NDROP  =  21 
DO  12  I  =  1 , NDROP 
GG  =  <1-1 >/3 , 0 

12  D I  AM  0<  I  >  =  5.5505428  *  <2.0**GG>  *  1  ,  OE-6 
C  STEP  SIZE  FOR  TARGET  VALUE  AND  FINAL  TARGET  VALUE  < M > 

DATA  XPT/. 00500/ 

DATA  XPTT/. 05635/ 

DATA  XPS/. 00135/ 

DATA  11/0/ 

DATA  KK/0/ 

C  MASS  FRACTION  OF  FUEL  IN  AMBIENT  AIR  SURROUND  INC  DROP 
DATA  YF 1/0.0/ 

C  FUEL  INLET  TEMPERATURE  < K ) 

DATA  FTEMPO/31 1 . / 

C  BOILING  POINT  OF  FUEL  <K> 

DATA  TBN  /560.0/ 

C  EXPANSION  COEFF  OF  FUEL 
DATA  CEX/ 0.00055/ 

C  CONE  ANGLE  OF  SPRAY  < DEGREES) 

DATA  PHI/45./ 

C  DISTANCE  FROM  NOZZLE  TO  SHEET  BREAKUP  <M> 

DATA  DBREAK/0. 005/ 

C  SWIRL  COMPONENT  OF  SPRAY  < DEGREES) 

DATA  FSWIRL/0 . 0/ 

C  INITIAL  SPRAY  VELOCITY  <M/3> 

DATA  VSPRAY/20 . 76/ 

C  MOLECULAR  WEIGHT  OF  AIR 
DATA  MWA/28,96/ 

C  INLET  AIR  VELOCITY,  X,Y,Z  COMP  < M/SEC) 

DATA  UA, VA, WA/08 ,54, 0 . , 0 . / 

C  STEP  SIZE  FOR  TIME  < SEC ) 

DATA  OELTAT/O. 00001 0/ 

C  LATENT  HEAT  OF  VAPORIZATION  AT  NORMAL  BOILING  POINT  <KJ/KG> 

DATA  LT6N/224 . / 

C  FUEL  DENSITY  <KG/M**3> 

DATA  RHOF 0/776 . 6/ 

C  MOLECULAR  WEIGHT  OF  FUEL 
DATA  MWF/226 . 43/ 

C  NUMBER  OF  CARBON  ATOMS  IN  FUEL  MOLECULE 
DATA  NOChRB/ 16/ 


<  USED  TO  CALCULATE  VAPOR  PRESSURE  > 
R  2 


AIR  PRESS<  ATM  >  = 


SURFACE  TENSION  OF  FUEL  < NT/M ) 

DATA  SIGMA/28 . OE-3/ 

CRITICAL  TEMP  OF  FUEL  <K) 

DATA  TCR/720 ■ 6/ 

CRITICAL  PRESSURE  OF  FUEL  <KPA> 

DATA  PCR/1421 ./ 

CRITICAL  DENSITY  OF  FUEL  <KG/M**3> 

DATA  RHGCR/239 . 0/ 

CRITICAL  COMPRESSIBILITY  OF  FUEL 
DATA  ZCR/0 . 225/ 

ACENTRIC  FACTOR  FOR  FUEL  <USED  TO  CALCULATE  V 
SEE  API  DATA  BOOK ,  CHAPTER  2 
DATA  ACENTR/0 . 74 1 8/ 

WRITE  < 1 ,999) 

)99  FORMAT  <"ATEMP0<K)  «,  AIR  PRESS<ATM>  =  " > 
READ  < 1 ,  *)  ATEMP0,P 
WRITE  < 1 ,990) 

)9 0  FORMAT  < "  XBAR  < MICRONS)*,  N  * ,  “ > 

READ  < 1 ,  ♦>  XBAR, N 
XBAR  *  XBAR  *  1 . OE-6 
WRITE  <  6 , 996  >  ATEMPQ,P 

)96  FORMAT  <"ATEMP0<K>  »M,F6.0,M  AIR  PRESS<  A 
CONVERT  PRESSURE  TO  KPA 
P  *  P  *  1 01 .325 
TAI*  ATENPO 

CALCULATE  QUANTITIES  WHICH  ARE  INDEP  OF  DROP  S 
DO  14  I  =  1 , NDROP 
14  TIMEv'I)  =  0.0 
SET  INITIAL  DROP  ARRAY  VALUES 

XPO  *  DBREAK  *  COS<PHI  *  3.1415/360.) 

YPO  *  DBREAK  *  SIN<PHI  *  3.1415/360.) 

ZPO  *  0. 

SPRAY  VELOCITY  COMPONENTS 

UPO  *  VSPRAY  *  CGS<  PH  I  *  3.1415/360.)  *  CO 

VPO  *  VSPRAY  *  SIN<  PHI  *  3.1415/360.)  *  CO 

WPO  =  VSPRAY  *  SIN< F SWIRL  *  2*3.1415/360.) 
DO  15  I  «  1, NDROP 
D<  I  )  »  DI AMO1!  I  > 

FTEMP< I  )  =  FTEMPO 


AIR  PRESS<  ATM  > 


INDEP  OF  DROP  SIZE 


“ , F6 . 2  ) 


3.1415/360.  ) 
3.1415/360. > 


3.1415/360.  ) 
3,1415/360. > 


COS<  FSWIRL*2+3 . 1415/360 
COS<FSWlRL*2+3. 1415/360 


UP<  I  ) 
VP<  I  ) 
WP<  I  ) 
XP<  I  > 
YP<  I  ) 
2P<  1  ) 
R<  I  )  » 
R<  0  >  > 


=  UPO 
*  VPO 
=  WPO 
■  XPO 
=  YPO 
«  2P0 
EXP<  - 
1  .  0 


< DI AM0< I )/XBAR >**N  > ) 


R<  NDROP+ 1 )  *  0.0 
DO  90  I  *  0, NDROP 
FF0<  I  )  *  R<  I  >  -  R<  1  +  1  ) 

HO  CONTINUE 

FF0<-1  )  *  0.  0 
OUTPUT  INITIAL  VALUES 
XBARMU  -  XBAR  *  1E4 

WRITE  <6,95)  XBARMU, N, FTEMP 0,TBN, PH I , VSPRAY, UPO, VPO, WPO 
?5  FORMAT  <F6. 1 ,F5.2,2F6.0,F6, 1 ,F7. 0,3F7, 1 > 

WRITE  <6,96)  ATEMPO.UA, VA,WA 
)6  FORMAT  <  4F1 2 . 1  ) 

WRITE  <6,97)  DELTAT ,P,XPQ, YPO, ZPO 
i?  FORMAT  •  FB.5,F5. 1 ,3FS.2) 


0113 
0120 
0121 
0 1  22 
0123 
0124 
0125 
0126 
0127 
0128 
0129 
0130 
0131 
0132 
0133 
0134 
0135 
0136 
0137 
01  33 
0139 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 
0156 
0157 
0158 
0159 
0160 
0161 
0162 
0163 
0164 
0165 
0166 
0167 
0168 
0169 
0170 
0171 
0172 
0173 
0174 
0175 
0176 
0177 
0178 


C  COMPUTE  TSST,  THE  STEADY  STATE  DROPLET  TEMPERATURE  <K>;  ALSO  LAMBST 
C  THE  STEADY  STATE  EVAPORATION  CONSTANT  <h**2/S> 

99  CONTINUE 
TSL  =  200. 

TSH  =  T  C  R  - 1 . 

IF  <  TSH . GT . TAI >  TSH  =  TAI-1 . 

DO  1 000  1  =  1 , 1 00 
IF  < I .EG. 1 >  GO  TO  1010 
IF  < I . EG . 2  >  GO  TO  f  020 
TS  *  <  TSL  +  TSH  >  /  2. 

GO  TO  1030 
1010  TS  =  TSL 

GO  TO  1030 
1 020  TS  *  TSH 
1030  CONTINUE 

TR  *  TS  +  <1./3.)  *  <  TAI  -  TS  > 

L  =  LTBN  *  <  <  TCR-TS  >/<  TCR-TBN  >  >**G . 38 
C  THIS  EXPRESSION  TO  BE  USED  ONLY  WHEN  DATA  FOR  FOLLOWING  EXPRESSION 
C  NOT  AVAILABLE 

PFS  =101  . 325*EXP<  2. 30259  *<-16, 188  +  3;=:.  1 24*TS/TBN-2  0 . 704*<  TSXTPN  >+ 
C  +4 . 768*<  TSXTBN  >**3  >  ) 

TREDUC  =  TSXTCR 

A  =  36/TREDUC  -  35.  -  TREDUC**6  +  96 . 73+ALGG 1 0< TREDUC > 

B  =  AL0G1 0< TREDUC >  -  0.0364*8 
C  =  ?*ALGG1 0< TREDUC)  -  0.!18*A 
TERM  1  =  C  -  1.192*B 
TERM2  =  4,93+6 

PREDUC  =  EXP<  2 . 30259  *  <  TERM i  +  ACENTR  *  TERM2  > > 

PFS  =  PREDUC  *  PCR 

PD IFF  -  P  -  PFS 

IF  <PDIFF.LT, 0.0)  PDIFF  =  0,0 

YFS  =  PFS+MWF  X  < PFS+HWF  +  PDIFF+MWA ) 

YFR  =  YFS  +  < 1  ,X3. >*<YFI-YFS> 

YAR  =  1 .  -  YFR 

IF  <  TR . GT . 600 ,  )  GO  TO  I  040 

CPA  =  4.184  *  <0.244383  -  4.20419E-5  *  TR  +  9.61128E-8  *  TR*  =  2 
8.  -  1.16393E-11  *  TR++3) 

GO  TO  1 050 

1040  CPA  =  4.184  *  <0.200331  +  7.71027E-5  *  TR  -  S.56726E-9  *  TR+  +  2  - 
8.  4 , 75772E-  1  2  *  TR++3 > 

1 050  CONTINUE 

CPV  =  <  0,363  +  0 . 000467+TR >  *  <5.  -  0.001 +RHOFO ) 

CPG  =  YAR  *  CPA  +  YFR  *  CPV 

LEFT  *  <  PXPFS  )  -  <MOF/MWA)*LA.CPG+<TAI-TS)>  -1. 

IF  < I , LT . 3 )  GO  TO  1 060 

IF  <ABS< LEFTJ.LT. 0. 0C5  >  GO  TO  1033 

IF  <  LEFT ,  LT . 0 . 0 )  GO  TO  103  0 

TSL  *  TS 

GO  TO  1000 

1 080  TSH  =  TS 

1060  IF  <1 .EG. 1 .AND.LEFT.LT. 0. 0>  GO  TO  1099 
IF  < I , EG. 2 . AND. LEFT. GT . 0 . 0 >  GO  TO  1099 
1 000  CONTINUE 
1 083  WRITE<  6, 1 08) > 

1081  FORMAT < "  TS/  TSL,  TSH,  LEFT,  P,  PFS,  L,  CPG  "  ) 

WR I TE<  6 , *  )  TS, TSL, TSH, LEFT, P, PFS, L, CPG 

1 070  TSST  =  TS 

WRITE< 1 , 1 075)  TSST, LEFT 
1075  FORMAT <  FI 0.2.F15.3) 

IF  <  TR . GT . 85 0 ,  >  GO  TO  1 1 00 


0179  KA  *  IE-4  *  <  1.1 3941E-3  +  1  .  01  7Q9E-3h.TR  -  5 . 825??E-7*Tfi**<.  + 

0180  &  2 . 46855E- 1  0h«TRh.*3> 

0181  GO  TO  1110 

0182  11  00  KA  «  1.E-4  *  <-.0999992  +  1  .  1 9  034E-3h.TR  -  5 . 29096E-?*TR**2  + 

0183  &  1 . 1 1208E-1 0*TR**3 > 

0184  1110  CONTINUE 

0185  E  3  2  .  -  0, 0772*< TR/TBN >*+2 

0186  KV  *  1.E-6  +  <13.2  -  0.0313  *<TBN-273.>>  *  < TR/273 . >**E 

0187  KG  *  YAR  *  KA  +  YFR  *  KV 

0188  WRI TE<  1,1111) 

0139  1111  FORMAT< "  KG, KA, KV, YAR , YFR, E , TR  "> 

0190  WRITE< 1,*)  KG, KA, KV, YAR, YFR, E, TR 

0191  BM  =  <  YFS  -  YFI>  /  <1  -  YFS  > 

0192  BT  *  CFG  *  <TAI  -  TSST  >  /  L 

0193  RHOFT  *  RHOFO  *  <  1  -  1  ,  8*CEXH"<  TS -288 , 6  >  -  <  0 . 09  0*<  TS-288 , 6  >*  + 

0194  8.  ✓  <  T CR  -  283 . 6 >**2 > > 

0195  LAM6ST  *  8  *  ALOG< 1 +BM  >  /  <  RHOFT  *  < CPG/KG  1  > 

0196  C  DENSITY  OF  AIR  <KG/M**3> 

0197  2AIR  *1.0 

0198  101  H  *  < 0. 0031 03/TR>  +  P  /  ZAIR 

0199  2AIRN  =  <  1  .A  1  ,-H>)  -  <  7349  .  /TR*+ 1  .  5  >  *  <  HA  1  . +H  )  > 

0200  IF  <  A6S<  <  2AIRN-2AIR )/2AIRN  > . LT . 0 , 0005  >  GOTO  102 

0201  2AIR  *  ZAIRN 

0202  GOTO  101 

0203  102  CONTINUE 

0204  RHOA  -  MWA  /  ZAIR  *  < PA  8 . 3 1 4*TR  >  ) 

0205  C  DENSITY  OF  FUEL  VAPOR  <KGAI**3> 

0206  CALL  VPRDN  <  TCR , PCR , MWF , TR , P , RHOV  > 

0207  C  VISCOSITY  OF  AIR  < KG7M*S )  SUTHERLAND  CORRELATION 
0208  MUA  *  1.425E-6  *  SQRT<  TR  >  A1  +  <100/TR>> 

0209  C  VISCOSITY  OF  THE  FUEL  VAPOR  AND  MIXTURE  WITH  AIR  <KG/M*S> 

0210  CALL  VPRVS<  TR , TCR , 2CR , MUA, MU V, MUG , MWA , MWF , YAR , YFR  ) 

0211  RHOG  -  1.  /  <  YAR/RHOA  +  YFR/RHOV > 

0212  WRITE< 1 , 11 12 > 

0213  1112  FORMAT< "  LAMBST , BM, BT , RHOFT , CPG , KG  "> 

0214  WRITE<1,*>  LAMBST, BM,BT, RHOFT, CPG, KG 

0215  C 
0216  C 
0217  C 

0218  C  COMPUTE  DROP  HEAT-UP  PARAMETERS 

0219  C  COMPUTE  QUANTITIES  WHICH  DEPEND  ON  DROP  SIZE 

0220  C 

0221  98  CONTINUE 

0222  C 

0223  DO  155  I*1,NDROP 

0224  IF  <TINE<  I  >.LT,  DEL  TAT  >  TSi><  I  >  =  FTEMPu 

0225  100  IF  <  <  TSD<  I  >-FT£MP0  >  .  GE  .  <  0 . 95h<<  TSST -F TEMP 0  >  >  >  GOTO  3060 

0226  TR  -  TSD< I >  +  <1./3.)  *  <TAI  -  TSD< I  )  > 

0227  C  DENSITY  OF  AIR  <KG/M**3> 

0228  ZAIR  -1.0 

0229  103  H  -  <  0. 0031  037TR)  *  P  /  ZAIR 

0230  ZAIRN  »  <t,A1.-H)>  -  <  7349 . 7TR*+  \  .  5  >  +  <  HA  1  .  +H  )  > 

0231  IF  < AB3<<ZAIRN-ZAIRVZAIRH>.LT. 0. 0005 >  GOTO  105 

0232  ZAIR  -  ZAIRN 

0233  GOTO  1 03 

0234  105  CONTINUE 

0235  RHOA  -  MWA  /  ZAIR  *  < PA  8 . 3 1 4  +  TR > > 

0236  C  DENSITY  OF  FUEL  VAPOR  <KG2M**3> 

0237  CALL  VPRDN  < TCR , PCR , MWF , TR , P , RHOV > 

0238  C.  VISCOSITY  OF  AIR  <KG/H*S»  SUTHERLAND  CORRELATION 


0239  MUA  -  1  . 425E-6  *  SQRKTR)  /<  1  +  <100/TR>> 

0240  C  VISCOSITY  OF  THE  FUEL  VAPOR  AND  MIXTURE  WITH  AIR  <KG'M*S> 

0241  CALL  VPRVS<TR,TCR,ZCR,MUA,MUV,MUG,MWA,MWF, YAR, YFR) 

0242  TREDUC  =  TSD< I  VTCR 

0243  A  =  36/TREDUC-  -  35 .  -  TREDUC* +6  +  96 . 73*AL0G 1 0< TREDUC > 

0244  B  =  ALOG1 0< TREDUC >  -  0.0364*A 

0245  C  =  7*AL0G 1 0<  TREDUC )  -  0.118+A 

0246  TERM1  =  C  -  1.192*B 

0247  TERM2  *  4 , 93+B 

0248  PREDUC  =  EXP-i  2.30259  *  <  TERM!  ♦  ACENTR  *  TERM2  >  > 

0249  PFS  -  PREDUC  *  PCR 

0250  PDIFF  =  P  -  PFS 

0251  IF  <PDIFF.LT. 0.0)  PDIFF  -  0.0 

0252  YFS  =  PFS*MWF  X  <PFS*MWF  +  PDIFF*MUA> 

0253  YFR  =  YFS  +  <1./3.)  *  <YFI-YFS> 

0254  YAR  ■  1 .  ~  YFR 

0255  IF  <  TR . GT . 600 .  )  GO  TO  3040 

0256  CPA  =  4.184  *  <0.244368  -  4.20419E-5  +  TR  +  9.61128E-8  *  TR*+2 

0257  *■  -  1  ,  1 6383E- 1  1  *  TR**3  > 

0258  GO  TO  3050 

0259  3040  CPA  =  4.184  *  <0.208831  +  7.7J027E-5  *  TR  -  8.56726E-9  *  TR**2  - 

0260  &  4 . 75772E- 1 2  *  TR**3 > 

0261  3050  CONTINUE 

0262  CPV  =  <0.363  +  0.000467+TR)  *  <  5 .  -  0.001*RHOF0> 

0263  C'PG  *  YAR  *  CPA  +  YFR  *  CPV 

0264  IF  <  TR , GT . 85 0 .  >  GO  TO  3100 

0265  KA  =  1.E-4  *  <1.13841E-3  +  1 . 0 1 739E-3*TR  -  5 . 82577E-7*TR**2  + 

0266  &  2. 46855E-1 0*TR**3> 

0267  GO  TO  3110 

0268  3100  KA  *  1.E-4  *  <-.0999992  +  1 . 19034E-3*TR  -  5 . 29 096E-7*TR**2  + 

0269  &  1 . 1 120CE-10*TR**3> 

0270  311 0  CONTINUE 

0271  E  »  2.-  0. 0372*<TR/'T6N>**2 

0272  KV  »  1.E-6  *  <13.2  -  0.0313  *< TBN-273 .  > >  +  < TR7273 . >**E 

0273  KG  »  YAR  *  KA  +  YFR  *  KV 

0274  L  =  LTBN  +  <  <  TCR-TSD< I ) >/< TCR-TBN  >  >**0 . 38 

0275  BT  *  CPG  *  <TAI  -  TSST )  /  L 

0276  BM  *  <  YFS  -  YF1>  /  <1  -  YFS) 

0277  RHOG  -  1.  /  < VAR/RHOA  +  YFR/RHOV  > 

0278  RHOFT  -  RHOFO*< 1 -1 . S*CEX*< TSD< I  >-288 . 6 )  -  <  0 . 090*< T3D< I >-263 . 6 >* 

0279  S.  X  <  TCP  -  288.6)**2>) 

0280  LAMBHU< I  )  -  3  *  ALOG< 1 +BM  )  /  < RHOFT  *  <CPG/KG>> 

0281  C  PRANDTL  NO. 

0282  3060  PRNDT  -  CPG  *  MUG  7  KG 

0283  C  REYNOLDS  NO, 

0284  RE<  I  )  -  D<  I  )  *  RHOG  SQRT<  <  UP<  1  >-UA  >**2  +  <  VP<  I  >  -VA  )**2  + 

0285  t  <WP<I)-WA>**2>  7  MUG 

0286  C  MODIFY  REYNOLDS  NO,  USED  TO  CALCULATE  MASS  LOSS  TO  INCLUDE  A 

0287  C  TERM  FOR  TURBULENCE  INTENSITY  EQUAL  TO  20  PERCENT  OF  X-AXIS  VEL 

0288  REM<  I  )  -  CXI)  *  RHOG  *< SQRT< < UP< I >-UA >**2  +  <  VP<  I  >-VA  >**2  * 

0289  &  <WP<  I  >-WA>**2>  ♦  0 . 2*UA  >  /  MUG 

0290  IF  <  <  TSD< I )-FTEMP0 ) . GE . <  0 . 95*<  TSST-FTEMPO  ) )  )  GOTO  3070 

0291  LAM8A< I >  -  LAMBHU< I  >  *  <1  +  0,276  *  REIKI  >**0.5  *  PRNDT** 0.33) 

0292  GOTO  3080 

0293  3070  CONTINUE 

0294  LAMBA<  I  )  »  LAMBST  *  <  1  +  0,270.  *  REIKI  >**0.5  *  PRNDT**  0 , 33  .> 

0295  3080  CONTINUE 

0296  MASSD  ■  <3.141596  X  6>  *  RHOFT  *  D< I >**3 

0297  MDOT  *  <3.141596  X  4 . >  *  RHOFT  *  LAM3A< I >  *  DU) 

0298  T  -  TSD< I >  /  TBN 


0299  CPF  -3 . 52390-6 . 15627*T+6 , 99895*T**2-2 . 4?850*T+*3+NOCaRB*< -0 . 016019 

0300  &  -0. 003128*T*0, 1?2136*T**2- 0. 058406*T**3> 

0301  TSD< I >-TSD< I >+<MDOT  *  L  /< CPF  *  MASSD >  ♦  <<6T/6M>  -  1.)  *  DELTA! > 

0302  IF  <  <  T8D< I >-FTEMP0  > . GE . <  0 . 95*<  TSST-FTEMPO  >  >  >  TSD< I >-TSST 

0303  D<  I  )  -  0<I>  -  <  LAMBA<  I  )/<  2  .  +&<  I  >  >  >  *  OELTAT 

0304  IF  <D<  I  >.LE.  0.9E-6)  D< 1 >-0 . 999E-6 

0305  C  DRAG  COEFFICIENT,  CD 

0306  IF  <RE< I >.GT .80. >  GO  TO  110 

0307  CO< I >  -  27.  *  RE< I >**< -0 . 84 j 

0308  GO  TO  130 

0309  1  1  0  IF  <RE<  I  ).GT  .  1  0000)  GO  TO  120 

0310  CD<I>  -  0.271  *  RE< I >**0.217 

03 i 1  GO  TO  130 

0312  120  CD<I>  =  2. 

0313  130  CONTINUE 

0314  C  DRAG  FACTOR,  F 

0315  F<  I  >  -  18.  +  MUG*CD< I ^  *  RE< I >  /(RHOFT  *  D< I >*♦2  *24.) 

0316  C  NEW  DROPLET  VELOCITY  AND  POSITION 

0317  UP< 1  )  *  -F< I >  *  <UP<I>  -  UA>  *  DELTAT  +  UP< I > 

0319  VP<  I  >  -  -F<  I  >  *  VP<I>  *  DELTAT  +  VP<  I  > 

0319  WP<  I  >  -  0. 

0320  XP< I  )  =  UP<I>  +  DELTAT  +  XP< I > 

0321  YP<  I  >  *  VP<I>  *  DELTAT  +  YP<  I  > 

0322  2P< I  >  -  0. 

0323  150  CONTINUE 

0324  IF  <D< I >.LE. 0.9E-6)  GO  TO  152 

0325  TIME<  I  >  =  TIMEU)  +  DELTAT 

0326  GO  TO  153 

0327  152  D<  I  >  «  0  9E-6 

0328  KK  =  I 

0329  GO  TO  155 

0330  153  CONTINUE 

0331  IF  <XP< I J.LT.XPS)  GO  TO  100 

0332  155  CONTINUE 

0333  C 

0334  C 

0335  C 

0336  DO  211  I  *  1 , NDROP 

0337  211  DMIK  I  )  -  D< I  >  *  1 . 0E6 

0333  GOTO  312 

0339  DO  312  1-1 .NDROP 

0340  WRITE  <6,310)  I ,  TIME<  I  >,  DMtK  I  >,  XP<  I  >,  FF0<  I  > 

0341  310  FORMAT  < IX, 13, FI 0,5,F8.2,F9.2,F9.6> 

0342  312  CONTINUE 

0343  C  COMPUTE  ROSIN-RAMMLER  PARAMETERS  OF  NEW  DISTRIBUTION 
0344  D<  KK  )  -  0.9E-6 

0345  D<  0)  ■  D< 1 > 

0346  DIAM0<  0  >  -  DIAM0< 1 > 

0347  UP<  0)  ■  UP<  1  ) 

0348  D<  NDROP+ 1 )  -  D< NDROP > 

0349  DIAMQ<  NDROP ♦ 1 >  =  DIAM0< NDROP > 

0350  UP<  NDROP* 1  >  -  UP< NDROP > 

0351  SUMFF  -  0,0 

0352  DO  160  I -KK, NDROP 

0353  160  S'JMFF  -  SUMFF  FF0<  I  )  *<<D<I>  /  DI  AM0<  I  >  ) +*3  /  UP<  I  )  +  (CKI+lJ 

0354  «,  /O I  AM0<  I  ♦  1  )  >+*3  /  UP<  !♦  1  >  >22 

0355  DO  170  I-KK, NDROP 

0356  170  FF<  I  )  -  FF O',  I  )  *  <  <  D<  I  >/D  I  AM0<  I  )  >  +  *3  2  UP<  I  ">  +  <  D<  I  + 1  VPIAM0<  I  ♦ 

0357  «.  1  )  )*"*3  /  UP<I+1>>  /  2  /SUMFF 

0358  DO  175  I-KK, NDROP 


035*  17S  R<I)  -  0. 

0360  DO  1  80  I  <*KK ,  NDROP 

0361  DO  130  J-I, NDROP 

0362  180  R<  I  )  ■  R<  I  >  ♦  FF<  J  ) 

0363  K  *  NDROP 

0364  DO  181  I=NBROP , KK+2, -1 

0365  181  IF  <R<  I  >.LT.  0.  OOt  >  K  ■  1-1 

0366  KL  =  KK  +  1 

036?  DO  190  1=KL,K 

0368  ZDATA< I )  «  ALOG< 1 XR< I >  ) 

0369  YDAT A<  I  >  =  ALGG< ZDATA< I  )  ) 

0370  190  XDATA<  I  >  =  ALOG<D<I>> 

0371  SUMX  -  0. 

0372  SUMY  -  0. 

0373  SUMX 2  *  0. 

0374  3UMY2  =  0. 

0375  SUMXY  *  0. 

0376  DO  200  I»KL,K 

0377  SUMX  *  SUMX  +  XDATA< I > 

0373  SUMY  «  SUMY  +  YDATA< I > 

0379  SUMXY  =  SUMXY  +  <  X0ATA< I )  *  YDATA< I > ) 

0380  SUMX2  *  SUMX2  +  XDATA< I  >**2 

0381  200  SUMY2  -  SUMY2  +  YDATA< I )**2 

0382  KD  =  K  -  KL  +  1 

0383  M  =  < SUMXY  -  < SUMX+SUMYXKO > >  X  < SUMX2  -  < SUMX**2XKD > > 

0384  B  =*  <  SUMY  -  M  +  SUMX  >/KD 

0385  SIGMAX  »  SQRT<<SUMX2  -  < SUMX**2 >XKD >  X  <KD-1>) 

0386  SIGMAN  -  S0RT<<SUMY2  -  < SUMY**2  >XKD >  X  <KD-1>) 

038?  CC  -  M  *  SIGMAX  X  SIGMAY 

0383  NNEW  »  M 

0389  XBARNW  -  EXP< -BXNNEW  ) 

039  0  XBRNWM  -  XBARNW  *  1  .  0E6 

0391  IF  < NNEW , LT . 1 . 0 >  GOTO  202 

0392  Y*2 . - 1 . XNNEW 

0393  G  -  Y*’*Y*EXP<  -Y  >*SGRT<  2  .  *3 . 141  5XY  >*<  1  ,  + 1  .  X<  1  2  .  *Y  >+ 1  ,  X<  288  .  *Y*+2  > 

0394  it  -139.X<51840.*Y**3>-5?1 .X<2488320.*Y+*4 >) 

0395  SMD  *  XBRNWM  *  <1.-1. /NNEW)  X  G 

0396  GOTO  204 

0397  2  02  WRITE*'  1,203) 

0398  203  FORMAT  < "  VALUE  OF  N  LESS  THAN  1.0  ") 

0399  204  CONTINUE 

0400  C  OUTPUT  STATEMENTS 
0401  WRITE  <6,205) 

0402  205  FORMAT  <//,'  R-R  XBAR  < UM >  R-R  N  CORR . COEF . ' 

04  03  it  '  SMD  ') 

0404  WRITE  <6,210)  XBRNWM , NNEW , CC , SMD  ■ 

04  05  21  0  FORMAT  < FI  5. 1 ,F 15. 2, FI  5. 3, FI  5. 2) 

0406  WRITE  <6,215) 

0407  215  FORMAT  <1X,X,7H  DIAM  ,8H  TIME  , 7H  X  DR0P,7H  Y  DROP , 7H  Z  DROP, 

0408  C  8H  U-VEL , 8H  V-VEL.8H  W-VEL , 6H  FTEMP , 6H  RE) 

0409  WRITE  <6,220)  <DMU<  I  >,  TIME<  I  >,  XP<  I  >,  YP<  I  ).  ZP<  I  >,  UP<  I  >,  VP<  I  >,  WP<  I  > 

0410  C  TSD<  I  ),  RE<  I  ),  1-1,  NDROP  ) 

041  1  220  FORMAT  <  F? . 1 ,F8 . 5 , 3F7 . 4, 3F8 , 3 , F6 . 1 , F / .  1  > 

0412  223  CONTINUE 

0413  GOTO  2402 

0414  WRITE  < 1 ,2400) 

0415  24  00  FORMAT< "  LAMBHIK  I ),  LAMBA< I >,  TSD< I  )  "  ) 

0416  WRITE  <  1 ,2401  )  <  LAMBHU< I  ), LAMBA< I >, TSD< I  ) ,  I » 1 , NDROP  ) 

0417  2401  FORMAT  < 2E 1 5 . 5 , FI 5 . 3  ) 

0418  2402  CONTINUE 


0419 

0420 

0421 

0422 

0423 

0424 

0425 

0426 

0427 

0423 

0429 

0430 

0431 

0432 

0433 

0434 

0435 

0436 

0437 

0438 

0439 

0440 

0441 

0442 

0443 

0444 

0445 

0446 

0447 

0443 

0449 

0450 

0451 

0452 

0453 

0454 

0455 

0456 

0457 

0458 

0459 

0460 

0461 

0462 

0463 

0464 

0465 

0466 

0467 

0468 

0469 

0470 

0471 

0472 


IF  <  XPS . L£ .XPTT  )  GO  TO  230 
GO  TO  240 

230  XPS  =  XPS  +  XPT 
GO  TO  98 

1  099  UR ITE<  1 ,  1 098 > 

1093  FORMAT< "  PROBLEM  IN  EVALUATING  TSST  AT  ORIG  LIMIT  VALUES  " > 

UR ITE< 1 ,*  >LEFT,P,PFS,MUF,MUA,L,CPG,TAI,TS 
240  CONTINUE 
STOP 
END 

SUBROUTINE  VPRDN  < TCR , PCR , MWF , TR , P , RHOV  ) 

REAL  MWF 
R  =  8.314 

A  *  0.4278  *  R**2  *  TCR**2.5  /  PCR 
L“  =  0.0867  *  R  *  TCR  7  PCR 
BB  -  &  /  <R*TR> 

AB  -  A  7  <B*R*TR+*t .5) 

2=1.0 

2210  H  =  BB  *  P  /  2 

ZN  *  < 1 . /< 1 . -H  >  >  -  AB*<  H/< 1 . +H  >  > 

IF  <  <  <  2N-2  )/ZN  > . LT . . 00  05  >  GO  TO  2200 
2  *  2N 
GOTO  2210 
2200  2  *  2N 

RHOV  -  <  MWF/2  >  *  <  PA  R+TR  >  > 

RETURN 

END 

SUBROUTINE  VPRVS  <  TR, TCR. ZCR , MUA , MUV, MUG, MWA, MWF , VAR , VFR  > 
DIMENSION  A<  7 ) 

REAL  MUA , MUV , MUG , MUA , MWF ,OMEG, OMGINV 
EK  -  65 . 3+TCR*ZCR**3 . 6 

SIG  -  0. 01866*<MUFARHOCR/1  000.  >>**0.333  /  2CR**1.2 

A< 1 >»9. 1426362E-1 

A<  2  >*-t . 068936  0 

A<  3  >=6 . 8077797E- 1 

A<  4  >=-2 . 1 208677E- 1 

A<  5  >*3 . 44871 86E-2 

A<6>— 2.8188225E-3 

A<  7  >=9 . 1 590342E-5 

OMGINV  »  0. 

DO  2300  1*1,7 

2300  OMGINV  -  OMGINV  +  A<  I  >*  <  ALOG 1  0<  Tft/EK  >  >**I 
OMEG  -  1 ./OMGINV 

C  VISCOSITY  OF  FUEL  VAPOR  <PA*SEC  OR  KG/< M*SEC > > 

MUV  *  2.669E-8  •*<  MWF*TR  >**0 , 5  /  <SIG**2  *  OMEG) 

C  CALCULATE  VISCOSITY  OF  MIXTURE 

PH1 1?**.  1  .  +<  NUh/MUV  )**G . 5*<  MWV/MUA  >**0 . 25  >**2/<  8  ,  *<  1  .  +MWA/MWV  )  >** 
«<  0.5 

PHI21  ■<  1  .  ♦< MUV/MUA  >**0 . 5*<  MWA/MWV  >**0 . 25  >*+?/<  8  .  *<  1  .  ♦MWV/MUA  /  >** 

t  0.5 

MUG=MUA/<  1  .  +PH1 1 2*<  YFR/MWF/<  YAR/MWA  >  >  >  +  MUV/-;  1  .  +FH 1 2 1  *<  YAR/MWA/ 
&  <  YFR/MWF  > ) ) 

RETURN 

END 
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